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ABSTRACT 

Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been 
detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio beams and more 
compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than 
non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected 
MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio 
profiles using a single- altitude hollow-cone beam, with a core component when indicated by polarimetry; however, 
for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase, we assume that 
the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence 
intervals are determined using a maximum likelihood technique. We divide the light curves into three model classes, 
with gamma-ray peaks trailing (Class I), aligned (Class II), or leading (Class III) the radio peaks. Outer gap and slot 
gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit 
with pair- starved polar cap models. Distinguishing between the model classes based on typical derived parameters 
is difficult. We explore the evolution of the magnetic inclination angle with period and spin-down power, finding 
possible correlations. While the presence of significant off-peak emission can often be used as a discriminator 
between outer gap and slot gap models, a hybrid model may be needed. 

Key words: acceleration of particles - gamma rays: stars - pulsars: general 
Online-only material: color figures, figure sets 


1. INTRODUCTION 

Millisecond pulsars (MSPs; first discovered by Backer et al. 
1982) are thought to be old pulsars, spun up to extremely 
short periods (P < 25 ms) via accretion from a companion 
(e.g., Alpar et al. 1982), and they are often dubbed “recycled” 
pulsars. The recycled pulsar scenario is supported by the fact 
that ~80% of MSPs are in binary systems and the detection 
of millisecond X-ray pulsations from neutron stars in low-mass 
X-ray binaries (LMXBs; e.g., Wijnands & van der Klis 1998; 
Chakrabarty 2005), presumed to be the progenitors of radio 
MSPs. Further evidence supporting this model was provided by 
the discovery of a radio MSP that had shown LMXB behavior, 
and no pulsations, in the past (Archibald et al. 2009), and it is 
thought to be a missing link in the LMXB-to-MSP evolutionary 
chain. More recent observations of PSR J1 824—24521 in the 
globular cluster M28 (NGC 6626) transitioning from rotation- 
powered radio pulsar to accretion-powered X-ray pulsar, and 
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back, have made this interpretation even more certain (Papitto 
etal. 2013a, 2013b). 

The rotation periods of MSPs are observed to be increasing 
at a much slower rate than non-recycled pulsars (P typically 
~10“^^ s s“^ for MSPs versus "^10“^^ s s“^ for non-recycled 
pulsars, Lorimer & Kramer 2004). This leads to weaker inferred 
surface magnetic fields (^surf = P P)^^^ /(In < 

10^ G, for an orthogonal rotator assuming dipole spin down and 
with 7 ?ns the neutron star radius, c the speed of light in vacuum, 
and I the neutron star moment of inertia). Lee et al. (2012) have 
empirically defined MSPs as those pulsars satisfying 



The light-cylinder radii of MSPs (cylindrical radius where 
co-rotation with the neutron star requires moving at the speed 
of light, Ri^c = o/Q with Q = 27r/P) are on the order of 
tens to hundreds of kilometers as opposed to many thousands of 
kilometers in non-recycled pulsars. The radio beams and polar 
cap sizes of MSPs are also very broad, making them detectable 
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over a larger range of viewing geometries than non-recycled 
pulsars. Therefore, MSPs provide excellent opportunities to 
study the global pulsar magnetosphere in detail through analysis 
of pulse profiles at different wavelengths. 

Using timing solutions from radio observatories around the 
world (Smith et al. 2008), MSPs have been established as a 
class of high-energy (HE; ^0.1 GeV) emitters (e.g., Abdo et al. 
2009a, 2010b, 2010d; Guillemot et al. 2012b) via observations 
with the Large Area Telescope (LAT; Atwood et al. 2009), the 
main instrument on the Fermi Gamma-Ray Space Telescope. 
Additionally, steady point-source emission has been detected 
from the vicinity of more than a dozen globular clusters (Abdo 
et al. 2010a, 2009b; Kong et al. 2010; Tam et al. 2011; Nolan 
et al. 2012) consistent with emission from the combination 
of many MSPs. HE pulsations have been detected from two 
extremely luminous MSPs in globular clusters (Ereire et al. 
2011; Wu et al. 2013; Johnson et al. 2013). The population of 
known radio MSPs in the Galactic field^^ has been increased 
by ~50% through follow-up searches of unassociated LAT 
sources with pulsar-like characteristics (e.g.. Ransom et al. 
2011; Cognard et al. 2011; Keith et al. 2011; Ray et al. 2012; Barr 
et al. 2013; Bhattacharyya et al. 2013), suggesting that MSPs are 
generally gamma-ray emitters. Of these new radio MSPs, over 
75% are in binaries and 11 are “black- widow” systems (with 
extremely low-mass companions thought to have been ablated 
by the pulsar wind; see Roberts 2011, for a review), further 
supporting the recycling scenario. 

The origin of HE pulsed emission remains an important 
question in gamma-ray pulsar physics. It is commonly accepted 
that the observed HE gamma rays are primarily the result of 
curvature radiation from electrons/positrons accelerated along 
the magnetic field lines by the rotationally induced electric 
field, but the recent detections of pulsations from the Crab 
pulsar at energies up to "^400 GeV by VERITAS (Aliu et al. 
2011) and MAGIC (Aleksic et al. 2011, 2012) suggest that 
either an additional component is necessary or a different 
process is at work (e.g., Lyutikov 2012). The exact location 
in the magnetosphere where the acceleration occurs is still 
uncertain. The light curves of gamma-ray pulsars detected 
with the LAT strongly suggest that emission occurs in the 
outer magnetosphere rather than near the polar caps, in narrow 
gaps bordering the closed-field-line boundary, but it is not yet 
clear how the acceleration or emission is distributed. Several 
authors have attempted to address these questions by generating 
simulated light curves using either geometric models (e.g., 
Cheng et al. 1986b; Venter et al. 2009; Bai & Spitkovsky 2010a; 
Watters & Romani 201 1) or full radiation models (e.g., Harding 
et al. 2008; Wang et al. 201 1). 

We have simulated gamma-ray and radio MSP light curves 
assuming the vacuum retarded-dipole (VRD) magnetic field 
geometry of Deutsch (1955). These simulations have been used 
to fit the observed light curves of the 40 gamma-ray MSPs from 
which significant pulsed signals have been detected with the 
LAT in three years of sky- survey operations as reported in the 
second LAT catalog of gamma-ray pulsars (2PC hereafter; Abdo 
et al. 2013), all of which satisfy the inequality in Equation (1). 

We fit the simulated light curves to the observed profiles 
and estimate uncertainties on the best-fit model parameters, 
using either one or two-dimensional likelihood profiles, and 
we discuss the implications of trends in the best-fit parameters. 
Similar to Venter et al. (2012), we define three model classes 
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as follows: MSPs with gamma-ray peaks trailing the radio 
peaks (by <0.5 in phase) are Class I, those with gamma-ray 
and radio peaks nearly aligned in phase (to within 0.1) are 
Class II, and those with gamma-ray peaks leading the radio 
peaks (by between 0.3 and 0.1 in phase, corresponding to radio 
lags between 0.7 and 0.9 in phase) are Class III. In Appendix A, 
we provide, for each MSP, the observed and best-fit light curves, 
summarize the observational characteristics, discuss how the 
models match the data, and compare our results to those from 
other methods when possible. Appendices B and C provide 
confidence contours and maps of simulated emission on the sky 
for selected MSPs, respectively. 

Because of the choice of radio phase zero in 2PC, PSRs 
J0034— 0534 and J1810— 1744 have radio lags of 0.866 and 
0.849, respectively, but they have wide double-peaked radio 
and gamma-ray light curves with similar morphology, and 
thus are considered Class II for our purposes. Additionally, 
PSRs J1744-1134 and 12214-^3000 have radio lags of 0.2 to 
0.3 in phase but are considered Class III in our studies. Eor 
PSR J1744— 1134, this is due to the use of the “h” method 
for shifting the profile, which puts phase zero at the weaker 
radio interpulse. Eor PSR 12214-^3000, phase zero is placed at 
the peak radio intensity, the “p” method, and the radio lag is 
referenced to the first gamma-ray peak appearing later in phase, 
but we note that the highest gamma-ray peak occurs just before 
the highest radio peak. 

Espinoza et al. (2013) separated the gamma-ray MSPs known 
at the time into three types based on the characteristics of their 
gamma-ray and radio light curves. Their A-type MSPs, those 
that have the main gamma-ray peak aligned with the main radio 
pulse, overlap exactly with our Class II MSPs with the exception 
of PSR J1 8 10-1- 1744, which was not known to be a gamma- 
ray emitter prior to 2PC. Their N-type and W-type MSPs are 
a mix of our Class I (3) and III (1) MSPs. Both N-type and 
W-type MSPs are defined as those without the main gamma-ray 
and radio peaks occurring at or near the same phase, with the 
radio profile of the former dominated by a single pulse and the 
latter consisting of wide peaks covering most of the pulse phase. 
Espinoza et al. (2013) hypothesized that the W-type MSPs were 
possibly aligned rotators, explaining the wide radio profiles and 
high duty cycles. 

2. PULSAR GEOMETRY 

The viewing geometry of a pulsar is defined by two angles, 
that between the spin axis and the magnetic dipole (o') and 
that between the spin axis and the observer’s line of sight (O- 
In studies of radio pulsars, the geometry is often characterized 
using a and P = ^ — o', with the impact parameter, describing 

the closest approach between the line of sight and the magnetic 
axis. 

Using geometric simulations, it is generally possible to 
reproduce almost any light curve shape for some combination 
of model parameters and viewing geometry, especially when the 
number of free parameters is large. As such, it is important to 
use all the multi-wavelength information available to restrict the 
parameter space and more strictly test the models. Eor reasons 
discussed below, strong constraints on viewing geometry exist 
for only a few MSPs and other considerations guide our 
modeling. 

Pulsar viewing geometries can be constrained from radio po- 
larimetry by fitting the observed position-angle swings with pre- 
dictions from the rotating- vector model (RVM; Radhakrishnan 
& Cooke 1969). However, the position-angle swings of MSPs 
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are typically not well fit by simple RVM predictions resulting 
in weak or no geometric constraints (e.g., Xilouris et al. 1998; 
Stairs et al. 1999; Ord et al. 2004; Yan et al. 2011). Relativistic 
and higher-altitude corrections can be applied to the RVM and 
lead to better constraints in some cases (e.g., Blaskiewicz et al. 
1991; Craig & Romani 2012). 

Polarization measurements can be used to guide the type 
of radio model used (i.e., a cone and/or core beam following 
Rankin 1983). The characteristic S-shaped or nearly fiat linear 
polarization position-angle swings are indicative of a hollow- 
cone beam viewed near or far from the magnetic dipole axis, 
respectively. Sense-reversing circular polarization has been 
interpreted as evidence for a core beam, and we take it as 
such here. Han et al. (1998) note that this is not restricted 
to core beams, but when observed in conal components, the 
sense of circular polarization is correlated with that of the 
linear polarization, which is not true for core components. 
Additionally, Beskin & Philippov (2012) have shown that the 
sense reversal can be due to propagation effects as the radiation 
travels through the pulsar magnetosphere. Most Class II MSPs 
are observed to have little or no polarized emission, which may 
be indicative of caustic emission (see Section 3 and Dyks et al. 
2004; Venter et al. 2012). 

In binary systems, the inclination angle of the orbit can 
be constrained if significant Shapiro delay is measured (e.g., 
Demorest et al. 2010). This angle should approximate ^ well, 
assuming that the spin and orbital axes are at least nearly aligned, 
such that i ^ providing a constraint against which to test 
our results. For “black-widow” and “redback” systems (Roberts 
2011; Ray et al. 2012), observing radio eclipses indicates that 
our line of sight is near the orbital plane, and therefore, the 
pulsar spin equator, suggesting ^ > 60°. 

For non-recycled pulsars with observed X-ray pulsar wind 
nebulae, it is possible to constrain ^ by modeling the shape 
of the torus (e.g., Ng & Romani 2008). However, MSPs do 
not typically power bright pulsar wind nebulae as, given their 
advanced ages, there is little material for the particle wind 
to interact with. If thermal X-ray pulsations are detected, 
models of the hot spot on the polar cap can be used to fit 
the light curve and estimate a, and other parameters (e.g., 
Bogdanov et al. 2007). However, the resulting fits are typically 
unconstrained. 

For MSPs with helium white dwarf companions, Tauris & 
Savonije (1999) have derived a relationship between the mass 
of the white dwarf (Mwd) and the period of the binary orbit. With 
this estimate of Mwd and the binary mass function of a pulsar, 
it is possible to estimate the inclination angle of the orbital 
plane with respect to our line sight if the mass of the pulsar is 
known (e.g.. Guillemot & Tauris 2014). Even if the mass of the 
pulsar is unknown, knowing the orbital period and assuming 
reasonable values for the pulsar mass provides constraints on 
the orbital inclination angle. Then, assuming that the orbit and 
spin axes have aligned over time, this provides a constraint on 

Guillemot & Tauris (2014) have shown that ^ constraints 
derived in this manner tend to agree well with those from fitting 
the gamma-ray and radio light curves of MSPs as described in 
this paper (see Section 5) and others (e.g., Johnson 2011). 

While full radiation models have been used to model particu- 
lar LAT pulsar light curves (e.g., Hirotani 2011; Du et al. 2011, 
2013), such analyses are time consuming and require tailoring 
of parameters such as P and P for each pulsar. For fitting large 
numbers of pulsar light curves, geometric models prove to be 
much more practical (e.g.. Venter et al. 2009; Romani & Watters 
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2010; Johnson 2011; Kalapotharakos et al. 2012a; Pierbattista 
et al. 2014). 

3. EMISSION MODELS 

Polar cap models (e.g., Daugherty & Harding 1996) assume 
that gamma rays come from near the stellar surface above the 
magnetic polar cap. In non-recycled pulsars, such models predict 
that the gamma-ray flux above a few GeV should be strongly 
attenuated due to magnetic pair creation. However, observations 
with the LAT strongly disfavor such models (Abdo et al. 2010d, 
2013) as the primary source of HE photons. The inferred dipolar 
magnetic field strengths of MSPs are not strong enough to lead 
to a significant attenuation of the gamma-ray flux. However, 
comparisons of the observed gamma-ray and radio light curves 
of MSPs do not agree with standard polar cap models, which 
predict broad peaks from wide emission regions. 

Outer-magnetospheric emission models assume that the HE 
emission is concentrated at high altitudes out to Rlc- The most 
commonly used models, for MSPs, are the outer gap (OG; e.g., 
Cheng et al. 1986a) model, slot gap (SG; Muslimov & Harding 
2003, 2004) or two-pole caustic (TPC; Dyks & Rudak 2003) 
model, and the pair- starved polar cap model (PSPC; Harding 
et al. 2005). 

In the OG and SG models, the magnetosphere is assumed to 
be filled with the charge density of Goldreich & Julian (1969) 
except in narrow vacuum gaps bordering the surface of last- 
closed field lines (those that close at the light cylinder) where 
particle acceleration is possible. In the OG model, the gap is 
bounded below by the null-charge surface (NCS; the geometric 
surface across which the charge density changes sign, defined by 
the condition - B = 0). In the SG model, the gap extends from 
the stellar surface out to the light cylinder. In this paper, we will 
consider the TPC model to be a geometric representation of the 
SG model. These models all predict relatively narrow gamma- 
ray light curve peaks that appear at later phase than features in 
the radio profiles, assuming low-altitude radio emission. 

In the TPC and OG models, the narrow bright peaks observed 
in gamma-ray pulsar light curves are the result of caustic 
emission (Morini 1983). This occurs when emission from 
different altitudes on the trailing field lines arrives closely spaced 
in phase due to relativistic aberration and the finite speed of light. 

In the PSPC model, the charged particle density is not 
sufficient to screen the accelerating electric field over the 
entire open-field-line region making it available for particle 
acceleration out to high altitudes. This model is only viable 
for MSPs and non-recycled pulsars with relatively low values 
of E (i.e., those pulsars below the “pair-creation death line,” 
Harding et al. 2002, 2005). PSPC models predict broad gamma- 
ray peaks that may appear at earlier phase than those in the 
radio. 

With the first gamma-ray MSP detections it was clear that 
acceleration was occurring in narrow gaps for most sources 
(Venter et al. 2009). Harding & Muslimov (201 la) showed that 
an offset of the magnetic dipole axis from the center of the polar 
cap could lower the death line and lead to narrow peaks in MSP 
light curves, depending on the amount of offset. 

A subclass of gamma-ray MSPs has emerged in which the 
gamma-ray and radio peaks are (nearly) aligned in phase (e.g., 
Abdo et al. 2010b; Guillemot et al. 2012b), a phenomenon only 
observed in one non-recycled gamma-ray pulsar (the Crab, e.g., 
Abdo et al. 2010c). For this subclass of MSPs, Abdo et al. 
(2010b) introduced the altitude-limited TPC and OG models 
(alTPC and alOG, respectively) in which the radio emission is 
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assumed to originate in regions that are significantly extended 
in altitude and co-located with the gamma-ray emission regions. 
These models were refined by Venter et al. (2012) who noted 
that the caustic nature of the radio emission has interesting im- 
plications for the expected polarization properties. In particular, 
for most geometries, the caustic mixing of emission from dif- 
ferent altitudes leads to the observed emission being largely 
depolarized (Dyks et al. 2004). 

In addition to the altitude-limited models. Venter et al. (2012) 
explored a low-altitude slot gap (laSG) model in which emission 
occurs within a few R^s of the polar cap, from pair cascades 
along the SG inner edge. The polarization pattern expected for 
the laSG model is basically that of the RVM (Dyks et al. 2004), 
producing high levels of linear polarization with the classic 
S-shaped position-angle swing. 

Qiao et al. (2004) and Qiao et al. (2007) proposed the 
annular gap model (later expounded on by Du et al. 2010, 201 1, 
2012, 2013) to explain gamma-ray and radio emission from 
millisecond and young, non-recycled pulsars. In this model, the 
open volume of a pulsar magnetosphere is separated into two 
regions by the critical field lines, those which cross the light 
cylinder at the point of intersection with the NCS. The region 
between the magnetic axis and the critical field lines is called the 
core gap while the region between the critical and last-closed 
field lines is called the annular gap. If the binding energy on the 
surface of the neutron star is strong enough, only one or the other 
acceleration regions can form (only the core gap if > 0 and 
only the annular gap if • ft < 0, where Q, is the spin angular 
momentum vector and fJL is the magnetic axis vector). For low 
surface binding energy, both acceleration regions can form. The 
size of the annular gap grows with a decreasing spin period and 
is thus thought to be most important for millisecond and young, 
non-recycled pulsars known to emit gamma rays. If the annular 
gap is large enough, pair-production is possible, leading to the 
production of secondaries that are accelerated from the stellar 
surface to the NCS (or even beyond, possibly out to the light 
cylinder) where gamma rays are produced. While we do not 
simulate MSP light curves using the annular gap model, we do 
discuss comparisons with results of light curve fitting using this 
model by other authors (see Section 7). 

Petri (2009) proposed the striped wind model to explain 
HE pulsar emission. In this model, the observed gamma rays 
are produced via inverse Compton scattering of the cosmic 
microwave background off an electron-positron pair wind 
outside the light cylinder. This model was able to successfully 
reproduce the phase-resolved spectra of the Geminga pulsar as 
observed by EGRET. When comparing the simulated gamma- 
ray and radio light curves from this model, assuming a polar 
cap model for the radio beam, Petri (2011) found that the basic 
profile characteristics could be reproduced, but the predicted 
phase lags were usually too large. 

4. LIGHT CURVE SIMULATIONS 

We assume that the VRD magnetic field geometry of Deutsch 
(1955) is valid in the rest frame of an inertial observer, to first 
order in t/Rlc- For a given spin period and a, we begin by 
finding the rim of the polar cap, defined to be the contour on the 
surface of the neutron star from which the last-closed field lines 
emerge. This is done iteratively by selecting a magnetic polar 
angle (O'), which is close to the polar cap opening angle 0pc 
~ (QRns /c)F^, as well as a particular magnetic azimuthal angle 
(00 and then integrating along the magnetic field line originating 
at (^0 00 using a fourth-order Runge-Kutta routine. If the line 
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closes inside the light cylinder, the code chooses a smaller value 
of O', and vice versa. This is repeated until a field line is found 
which closes near the light cylinder, and this footpoint position 
is stored. The code then moves to the next (p', looking for the 
0' of the next last-closed field line, and continues this process 
until the rim of the polar cap has been fully defined. 

Since the VRD polar cap is not symmetric about the magnetic 
axis, a new coordinate definition is used to label positions on 
the stellar surface. A newer implementation of the so-called 
open- volume coordinates (rove, ^ovc) was introduced by Dyks 
et al. (2004), following initial work by Yadigaroglu (1997); 
Cheng et al. (2000). The “radial” coordinate is formally defined 
as Tovc = 1 ± doYc, where Jove is the minimum distance of 
a point from the polar cap rim, normalized by the polar cap 
radius Rpc = ^ns^pc ~ i^ similar to a 

scaled co-latitude which ranges from 0 to 1, labeling points (or 
cones) from the magnetic axis to the rim. However, given the 
asymmetric polar cap shape, contours of constant rove now label 
concentric, self-similar, deformed rings interior or exterior to the 
polar cap instead of circles of constant co-latitude on the stellar 
surface (see Figure 2 of Dyks et al. 2004). On the rim, rove = 1- 
The second coordinate /ove is analogous to the azimuthal angle, 
and measures arc length along a deformed ring of fixed rove 
in the direction of increasing azimuthal angle, with /ove = 0 
coinciding with zero azimuth. 

Rings of constant rove are defined between specified values 
of r^^“ and depending on the model and waveband to be 
simulated, as described in Sections 4.1 and 4.2. Magnetic field 
line footpoints are positioned equidistantly along the arc length 
of each ring (non-uniformly in azimuth). Therefore, there will 
automatically be fewer magnetic field lines positioned on the 
inner rings, since their circumferences are smaller, leading to a 
uniform placement of field lines over the PC such that one would 
not overestimate the emission originating at small colatitudes. 

We next follow a hypothetical electron /positron as it is 
accelerated along each field line between and r^^^ and 
collect the resulting emission in bins of pulse phase and 
The emission is assumed to be tangent to the local magnetic 
field line direction at the point of emission in a frame that co- 
rotates with the star. To calculate the phase and ^ at which 
an emitted photon would be observed, this local direction is 
first transformed from the inertial observer’s frame to the co- 
rotating frame (for details of this calculation see Johnson 2011) 
and used to calculate the emission direction as advocated by 
Bai & Spitkovsky (2010b). This direction is then transformed 
back to the lab frame correcting for relativistic aberration and 
time-of-fiight delays. 

The corresponding phase and ^ bins are incremented by a 
number of photons for gamma-ray models, either proportional 
to physically motivated emissivity profiles (see, e.g.. Venter et al. 
2009, 2012) or to the step length along the magnetic field line for 
uniform emissivity, and flux level for radio models (and divided 
by the solid angle spanned by each bin). Our models assume 
uniform emissivity along the field lines in the co-rotating frame 
except for the PSPC and laSG models. 

Simulated light curves are constructed by plotting the bin 
contents for a given ^ from the skymap generated for a given 
a. The ^ bins are defined with integer boundaries (e.g., [0°,1°), 
[1°,2°), etc.), but only the lower boundary is reported as the 
best-fit ^ in Section 7. Venter et al. (2009) have shown how the 
emission skymaps of models used to fit Class I and III MSPs 
change for different parameters, and Venter et al. (2012) have 
done the same for models used to fit Class II MSPs. Skymaps 
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corresponding to the best-fit geometries of select MSPs can be 
found in Appendix C. 

We assume a neutron star mass and radius of 1.4 Mq and 
7 ?ns = 10 km, respectively, which lead to / = 10^^ g cm^ 
assuming a spherical star with uniform mass density. These 
values only affect the PSPC model where we use an ana- 
lytic form of the accelerating electric field to calculate the 
number of emitted photons at each spatial step (Venter et al. 
2009). Different assumptions for the neutron star equation of 
state will mainly influence the stellar compactness (a 
Muslimov & Harding 2003), which in turn would change the 
electric field, the latter being roughly linearly dependent on the 
compactness. 

4.1. Gamma-ray Simulations 

For OG and TPC models, we assume uniform emissivity 
along the field lines in the co-rotating frame such that the number 
of emitted photons is proportional to the step length along the 
field line. For all gamma-ray emission models, particles are 
followed out to a radial distance of 1.2 7 ?lc (except for the 
altitude-limited models where this becomes a free parameter) 
but not beyond a cylindrical distance of 0.95 7 ?lc as the 
structure of the magnetic field near the light cylinder is not well 
known. 

We use a two-layer OG model in which particles are accel- 
erated in a vacuum gap between the surface of last-closed field 
lines (at the rim with rove = 1) and a specified r^^^ ^ 1. The 
particles are then assumed to emit HE photons via curvature 
radiation above the NCS on those field lines emerging from 
the stellar surface between r^^^ and a specified ^ in 
the standard OG model Following Wang et al. 

(2010), this emission layer should be small compared to the 
vacuum gap width (ujacc = 1 — and thus we constrain 

the size of this region (Wem) to be no more than half the size 
of the vacuum gap (i.e., ^ 0.5(1 - r^^^)). 

In the TPC model the emission layer and accelerating gap 
are the same. Particles are followed along field lines that 
originate on the stellar surface between r^^^ = 1 and some 
specified r^^^ ^ 1. For TPC models, the gap width is defined as 

1 _ mm 
ovc 

In PSPC models, the entire open volume is available for 
particle acceleration and emission is collected from all field 
lines with 0 ^ rove ^ 1- The accelerating field is expected 
to change drastically over this range of field lines, and thus 
the uniform emissivity assumption is no longer reasonable. 
Therefore, we use the solution of the accelerating field and 
emission prescription outlined in Venter et al. (2009) to calculate 
the number of photons emitted per step length along a field line. 

Class II MSP light curves are simulated using alTPC, alOG, 
and laSG models (Venter et al. 2012). The altitude-limited 
simulations are carried out similar to the standard OG and TPC 
models, with the exception that the maximum radial distance 
(^max) to which the emission is followed is a specified parameter, 
taking on values from 0.7 to 1.2 7 ?lc- 

In the laSG model, the emission peaks at a distance Sf above 
the stellar surface (typically 1 to 2 /?ns. Muslimov & Harding 
2003; Venter et al. 2012) and falls off exponentially above and 
below Sf with a characteristic length crin toward and cTout away 
from the star, both in units of R^s • We use this emission profile 
to modify the number of photons that would have been added 
to the skymap under the assumption of uniform emissivity. This 
emission profile is motivated by physical models involving 
particle acceleration, curvature emission, and pair cascade 
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formation. The acceleration and emission is taken to occur 
between = 0.95 0pc and a specified 

4.2. Radio Simulations 

A single-altitude hollow-cone and/or core beam is used to 
model the radio profiles of Class I and III MSPs. We follow 
the prescription of Story et al. (2007, who built on the work 
of Arzoumanian et al. 2002; Gonthier et al. 2004) but do not 
assume the same core-to-cone flux ratio. In particular, the two 
components are generated separately, and if both are used, the 
normalizations are fit separately. 

The polar cap is divided into rings, as described previously, 
between =1.2 0pc and =0.1 0pc (cone) or 0.0 (core). 
These values are chosen to allow the Gaussian tails of the beam 
components to be explored. The field line emerging from each 
ring segment is followed out to the emission height calculated 
by Kijak & Gil (2003): 



in units of 7 ?ns where v is the emitted frequency. Once this height 
is reached along a field line the observed phase and ^ values are 
calculated and that bin is incremented by the appropriate flux 
level (given in Harding et al. 2008). 

Equation (2) depends most strongly on P and v, but for typical 
values, the expected emission heights are <0.3 R^c Such mod- 
els do not produce radio peaks at the same phase as those from 
outer-magnetospheric emission models, and thus cannot explain 
the Class II MSPs. The flux of the core component is maximum 
along the magnetic axis, dropping off as a Gaussian with charac- 
teristic width (Tcore = l°5/(P/ls)^/^. The characteristic size of 
the cone beam is given by Pcone = H24(rKG)^'^^(T’/ls)“^/^ (the 
width at 0. 1% of the peak intensity. Story et al. 2007), which has 
overall and dependencies, consistent with Kramer 

et al. (1998) and Kramer et al. (1999), respectively. The flux 
of the cone component is a maximum at magnetic polar angle 
0^ = 0.52pcone and falls off as a Gaussian on either side with 
characteristic width cjcone = 0.18pcone (Harding et al. 2008). 

The radio profiles of Class II MSPs are simulated with alTPC, 
alOG, and laSG models. For the alTPC and alOG models, 
this implies that the radio emission can extend to substantially 
larger altitudes than given by Equation (2). In these models, 
the minimum radius of emission is also a parameter of the 
simulations. In particular, for alTPC models emission is only 
collected between a specified and alOG models 

the minimum radius of emission is taken to be max{/?^-j^, T^ncs}. 
where 7 ?ncs is the radius of the NCS and is a function of magnetic 
azimuth and a. 

4.3. Simulation Parameters 

We have generated simulations using the models described 
previously with parameter resolutions and ranges given in 
Table 1 for the TPC and OG models; Table 2 for the alTPC, 
alOG, and laSG models; and Table 3 for the PSPC model. When 
fitting light curves of Class II MSPs, the radio and gamma- 
ray emitting regions are allowed to have different gap widths. 
When using laSG models to fit the Class II MSP light curves, we 
require that (Tout > and require that all parameters, except the 
gap widths, be the same for the radio and gamma-ray emitting 
regions. 
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Table 1 

Simulation Parameter Resolutions and Ranges for Class I MSP Models 


Model 

P 

(ms) 

(X 

(o) 

(o) 

^acc 

(% 0Pc) 

^em 

(% 0pc) 

Resolutions: 

TPC 

1.0 

1 

1 

2.5 


OG 

1.0 

1 

1 

2.5 

2.5 

Ranges: 

TPC 

[1.5, 5.5] 

[1,90] 

[0, 180] 

[0.0, 10.0] 


OG 

[1.5,5.51 

[1,90] 

[0, 180] 

[0.0, 10.0] 

[0.0, 5.0] 


5. LIGHT CURVE FITTING 

We fit the gamma-ray and radio light curves by scanning 
over our model-parameter grid. We define phase zero in the 
simulated light curves to be the closest approach of the magnetic 
axis. When scanning over the parameter space, we introduce a 
phase shift (O, in units of light curve bins) of the simulated 
light curves as an additional fit parameter. For a given set of 
model parameters, we scan through all possible values of O 
(from OiOHy — 1 , with Uy the number of bins used in a given 
gamma-ray light curve) optimizing the normalizations of the 
simulated radio and gamma-ray light curves, on top of estimated 
background levels, separately. The gamma-ray light curves are 
fit with Poisson likelihood while the radio profiles are fit with a 
statistic and then the two statistics are combined as follows. 

In maximum likelihood techniques, problems are typically 
recast as minimizations of — ln(£), where C is the likelihood. 
For the gamma-ray light curves, we minimize 


-\n{Cy) = - In 


Uy—l 

n 

/=0 


4" ^y) ^^p{ 4" ^]/)} 

dyj ! 


(3) 


In Equation (3), dyj is the value of the ith gamma-ray light curve 
bin, by is the background estimate, and CyjX^ gives the value 
of the ith model light curve bin A/. The Cyj values are defined 
such that, for some reference bin 7 ^ 0, Cyj = "ki/Xf. When 
minimizing Equation (3), the only parameter that is optimized 
is . 

For the radio light curves we minimize 


0 5 

- 1 ii(£r) = — y]((cR,i/?>p + b^- jR,i)^ (4) 
;=o 


Table 3 

Simulation Parameter Resolutions and Ranges for Class III MSP Models 


Model 

P 

(ms) 

a 

(o) 

(o) 

Resolutions: 

PSPC 

1.0 

1 

1 

Ranges: 

PSPC 

[1.5, 5.5] 

[1,90] 

[0, 180] 


In Equation (4), Jr i is the value of the ith radio light curve bin, 
Z?R is the background estimate, ctr is the error used for each radio 
bin (see the following discussion), and gives the value 

of the ith model light curve bin Ri. The cr^ values are defined 
such that, for some reference bin R^r ^ 0, cr i = Ri/R^t. When 
minimizing Equation (4), the only parameter that is optimized 
is Rxp. We combine the gamma-ray and radio log-likelihood 
functions by adding Equations (3) and (4). 

When attempting a joint fit of the radio and gamma-ray light 
curves one difficulty arises due to the fact that the statistical 
uncertainty of the radio data is much less than that of the gamma 
rays. Without adjusting the radio uncertainty, this results in the 
likelihood essentially ignoring the information in the gamma- 
ray light curve when choosing the best-fit geometry. Given the 
simplistic nature of our radio model for Class I and III MSPs, 
this leads to inadequate solutions. 

To balance the relative contributions of the radio and gamma- 
ray profiles, we define an on-peak interval for the gamma-ray 
light curve and calculate the average relative uncertainty in this 
interval (cr^,ave)- We then calculate ctr = (r^ax x where 

^max is the value of the highest radio light curve bin, and use 
this as the absolute uncertainty for each radio bin. For those 
MSPs where we use more radio bins than gamma-ray bins, this 
uncertainty is decreased by the ratio of gamma-ray to radio bin 
numbers. 

The bins in the on-peak interval of the gamma-ray light curve 
are the dominant contribution to that part of the likelihood, and 
using (Tr as defined leads to a comparable contribution from the 
radio profile. However, this ctr is still somewhat arbitrary, so we 
have attempted to estimate systematic biases in our fits due to 
this choice (see Section 5.1). 

We construct likelihood profiles to estimate uncertainties on 
the best-fit model parameters. For a and ^ we use a two- 
dimensional likelihood profile to produce confidence contours 
that can be compared to those from polarization fits, otherwise 


Table 2 

Simulation Parameter Resolutions and Ranges for Class II MSP Models 


Model 

a 

( o ) 

( o ) 

^acc 

(% 0Pc) 

'^em 

(% 0pc) 

^max 

(Rlc) 

dR 

^min 

(Rlc) 

dR 

^max 

(Rlc) 

(Rns) 

<^out 

(Rns) 

Sf 

(Rns) 

Resolutions: 

alTPC 

1 

1 

2.5 


0.05 

0.05 

0.05 




alOG 

1 

1 

2.5 

2.5 

0.05 

0.05 

0.05 




laSG 

1 

1 

2.5 





0.1 

0.3 

0.2 

Ranges: 

alTPC 

[1,90] 

[0, 180] 

[0.0, 10.0] 


[0.7, 1.2] 

[0.14, 1.15] 

[0.2, 1.2] 




alOG 

[1,90] 

[0, 180] 

[0.0, 10.0] 

[0.0, 5.0] 

[0.7, 1.2] 

[0.14, 1.15] 

[0.2, 1.2] 




laSG 

[1,90] 

[0, 180] 

[0.0, 10.0] 





[0.1, 1.0] 

[0.3, 2.1] 

[1.2, 2.0] 


Notes. All alTPC, alOG, and laSG simulations are done using a 1.5 ms spin period. For the Wacc and Wem parameters the radio and gamma-ray 
simulations have the same resolutions and ranges. The lower limit for of 0.14 Rye corresponds to Rns for a 1.5 ms spin period. 
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we use one-dimensional profiles to estimate the uncertainties. 
In order to be conservative, we have chosen to report 95% 
confidence-level uncertainties. 

Our geometric models are simple when compared to the 
complex issue we are trying to address and the likelihood 
surfaces can be steep near the best-fit values leading to un- 
reasonably small uncertainties. Thus, when estimating the 
uncertainties, we multiply the log-likelihood differences by 
(^dof/2)/(— ln(£niax)), where £max is the maximum likelihood 
value and ^dof is the degrees of freedom in a particular fit. As- 
suming the log-likelihood differences follow a distribution 
(—A ln(£) = Ax^/2), this results in the best fit corresponding to 
a reduced x^ of 1. We use these rescaled log-likelihood differ- 
ences to estimate uncertainties on the best-fit model parameters. 

The confidence contours are often not simple shapes with 
confidence regions that are not simply connected. Thus, while 
we do report zb error bars for these angles in Section 7, when 
comparing to other constraints (such as RVM fits), it is best to 
use the confidence contours directly, provided in the auxiliary 
online material. 

5.7. Evaluation of Systematic Biases 

We have investigated multiple sources of systematic biases 
in our fitting procedure. The first is related to our choice of 
(Tr. Choosing a smaller (larger) value makes matching the radio 
profile more (less) important to the overall likelihood value and 
can thus affect the best-fit geometry. To assess the importance 
of this systematic, we refit a sub- sample of MSPs, one-third of 
each model class, varying the radio uncertainty by a factor of 
two. 

For Class I MSPs, changes in a and ^ , due to varying (Tr by a 
factor of two, were generally <30° and typically in different 
directions (i.e., if a increases ^ usually decreases and vice 
versa). There were a few extreme cases where the change in one 
parameter was as much as ~60°. In such cases, the best fit was 
jumping from one local maximum to another, but such changes 
are encompassed by the confidence contours (see Appendix B . 1 , 
for examples). Depending on the gamma-ray light curve, only 
fits with one of the models (TPC or OG) may be affected, but 
the general behavior of the two models was similar. The gap 
width parameters were generally unaffected, with any changes 
within the estimated statistical uncertainties. 

For Class II MSPs, varying ctr by a factor of two led to 
changes in ckf or ^ ^2°. The additional model parameters were 
either unchanged or changed by less than the estimated statistical 
uncertainties. 

For Class III MSPs, a and ^ typically changed by only a few 
degrees when varying ctr by a factor of two. However, for fainter 
MSPs in this class, the geometry was observed to change by as 
much as ~50°. 

The estimated background level in the gamma-ray light curves 
can also strongly affect the best-fit parameters. These values 
come from the weights estimated using spectral fits that can have 
sizable statistical uncertainty. From investigations in 2PC, we 
estimate that these background levels could be wrong by "^5%. 
This could change which model is preferred as it alters the off- 
peak emission level for which OG and TPC models give different 
predictions. To assess the importance of this systematic, we refit 
a selection of MSPs varying the gamma-ray background level 
by di5%. 
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For Class I and II MSPs, a modified background level led to 
changes in a and ^ typically less than a few degrees and ~10° 
at most. The most significant effect was in the — ln(£) values, 
which were seen to change by <7 in most cases. In one extreme 
case, for the Class II MSP PSR J 1939-1-2 134, which has a very 
high gamma-ray background level; the changes in the — ln(£) 
values were on the order of "^20. To be conservative, we only 
say one model is significantly preferred over another for a given 
MSP if the log-likelihood difference is ^15. Our models are 
not nested, so assigning a significance to any log-likelihood 
difference is not straightforward; however, we note that a log- 
likelihood difference of 15 amounts to the preferred model being 
more likely by a factor of e^^ (compare this to a 5a result for 
nested models with one degree of freedom corresponding to 
the preferred model more likely by a factor of ~e^^'^). As with 
changes in ctr, the gamma-ray gap width parameters were not 
significantly affected by this systematic. The emission altitude 
parameters in the fits of Class II MSPs were also not strongly 
affected, with any differences within the estimated statistical 
uncertainties. 

The best-fit geometries of Class III MSPs are insensitive to 
this systematic. If any change was observed it was ^1°. 

The VRD magnetic field geometry is only an approxima- 
tion. Charges will be pulled from the surface of the neutron 
star, populate the magnetosphere, and modify the field struc- 
ture. Attempts to construct an analytic model of a pulsar mag- 
netosphere, accounting for the presence of charges, have only 
been successful for aligned rotators (e.g., Michel 1974; Mestel 
& Pryce 1992). However, using magneto-hydrodynamic sim- 
ulations several authors have successfully produced numeric 
models assuming force-free conditions (i.e., no acceleration, 
Spitkovsky 2006) and for finite-conductivity spanning the range 
between vacuum and force-free (Kalapotharakos et al. 2012b; 
Li et al. 2012). One result of including charges is an increase 
in the size of the polar cap. This leads to larger predicted phase 
lags between the gamma-ray and radio pulse profiles (by as 
much as 0.1 in phase, Harding et al. 2011; Kalapotharakos et al. 
2012a). 

In order to evaluate how the use of the VRD approximation 
may skew the best-fit geometries, we repeated the fits for a subset 
of Class I MSPs allowing phase zero of the gamma-ray models 
to be different than the radio models by ^0.1. Additionally, we 
refit some of the Class III MSPs with OG and TPC models with 
this same strategy in an attempt to address if an increased phase 
lag allows the Class I models to explain the Class III light curves 
or if there really is evidence for pair-starved magnetospheres in 
these MSPs. 

For Class I MSPs, allowing the gamma-to-radio phase lag 
to increase beyond the VRD prediction led to changes in a and 
^ ^ 10°, almost always in the same direction (i.e. , both increased 
or decreased). In the case of one faint MSP, a was observed to 
change by ~30°. 

For Class III MSPs, allowing the gamma-to-radio phase 
lag to increase beyond the VRD prediction and fitting with 
TPC and OG models did not result in acceptable fits. The 
maximum allowed change of 0.1 was not sufficient to match 
the observations, suggesting that the light curves of these 
MSPs cannot be explained only by deficiencies in the VRD 
approximation. 

6. SOURCES AND DATA PREPARATION 

For each of the MSPs in Table 4, we have used the same 
LAT data set and radio profiles as 2PC. The P and spin-down 
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Table 4 

Parameters of 2PC MSPs 


PSR 

^obs 

Albs 

^obs 

Observatory^ 

E 

Blc 

d 

Ly/fa 

nylfo. 


(ms) 

(10“^^ s s“^) 

(MHz) 


(10^^ erg s“^) 

(10‘‘ G) 

(pc) 

(10^^ erg s“') 

(%) 

J0023+0923 

3.05 

1.08 

2000 

GET 

151.00 

5.99 

690!.2W 

4.56 ± 0.70 

3.00 ± 0.50 

J0030+045F 

4.87 

0.11 

1400 

NRT 

36.40 

5.82 

280+‘“ 

5.75 ±0.17 

15.8 ±0.50 

J0034-0534^ 

1.88 

0.29 

1400 

NRT 

173.00 

10.35 

540^‘i‘oo 

5.67 ± 0.41 

3.30 ± 0.20 

JOlOl-6422^ 

2.57 

0.44 

1400 

PKS 

101.00 

5.83 

550-Z 

3.79 ± 0.33 

3.80 ± 0.30 

J0102+4839 

2.96 

1.17 

1500 

GET 

175.00 

6.68 

oai q+500 

85.1 ± 10.2 

48.6 ±5.8 

J0218+4232^ 

2.32 

7.69 

1400 

NRT 

2432.00 

31.51 

2640'''^^^^ 

380 ± 20 

15.6 ±0.8 

J0340+4130 

3.30 

0.59 

1500 

GBT 

78.70 

3.62 

^ '■^'^-300 

72.9 ± 5.5 

92.6 ± 7.0 

J0437-4715'' 

5.76 

1.41 

1400 

PKS 

29.05 

1.39 

156 ± 1 

0.49 ± 0.03 

1.70 ±0.10 

J0610-2100^ 

3.86 

0.12 

1400 

NRT 

8.00 

1.10 

21540+5460 

98.5 ± 16.0 

1230 ± 201 

J0613-0200^ 

3.06 

0.87 

1400 

NRT 

120.00 

5.30 

900!** 

29.0 ± 1.8 

24.1 ± 1.5 

J0614-3329 

3.15 

1.78 

1400 

NRT 

220.00 

7.06 

1900-fo 

472 ± 12 

215 ±5 

J0751+1807^ 

3.48 

0.77 

1400 

NRT 

72.14 

3.62 

400!.™ 

2.54 ± 0.23 

3.50 ± 0.30 

41024-0719'“ 

5.16 

0.16 

1400 

NRT 

4.00 

0.62 

386!fg 

0.57 ±0.18 

14.2 ± 4.5 

J1 124-3653 

2.41 

0.58 

820 

GBT 

171.00 

7.83 


42.9 ± 4.7 

25.1 ±2.7 

J1 125-5825 

3.10 

6.09 

1400 

NRT 

805.00 

13.59 

zoiv_3YQ 

73.1 ± 16.4 

9.10 ±2.00 

J1231-141F 

3.68 

0.65 

1400 

NRT 

51.46 

2.89 

438!« 

23.6 ± 0.6 

45.9 ± 1.2 

J1446-4701 

2.19 

0.99 

1400 

PKS 

368.00 

13.03 

1460!™ 

19.0 ±3.7 

5.10 ± 1.00 

J1514-4946 

3.59 

1.87 

1400 

PKS 

160.00 

5.22 

940!™ 

47.6 ± 3.0 

29.7 ± 1.8 

J1600-3053^ 

3.60 

0.86 

1400 

NRT 

72.98 

3.51 

1630!2*,q 

16.9 ±8.9 

23.1 ± 12.2 

J1614-2230^ 

3.15 

0.30 

1400 

NRT 

37.79 

2.90 

650 ± 50 

12.3 ± 1.0 

32.6 ± 2.6 

J1658-5324 

2.43 

1.10 

1400 

PKS 

302.00 

10.61 

930+11^ 

29.9 ± 2.4 

9.90 ± 0.80 

J1713+0747^ 

4.57 

0.83 

1400 

NRT 

34.20 

1.90 

1050+“ 

13.4 ± 1.9 

39.2 ± 5.6 

J1741+135F 

3.75 

2.91 

1400 

AO 

217.60 

5.83 

1080!“ 

3.36 ± 1.13 

1.50 ±0.50 

J1744-1134^ 

4.07 

0.70 

1400 

NRT 

41.08 

2.33 

417!\*7 

6.76 ± 0.52 

16.5 ± 1.3 

J1747-4036 

1.64 

1.33 

1400 

PKS 

1160.00 

31.19 

3390!™ 

136 ±31 

11.7 ±2.7 

J1810+1744 

1.66 

0.46 

2000 

GBT 

397.00 

17.85 

2000!™ 

112± 8 

28.2 ±2.1 

J1823-3021A 

5.44 

338.00 

1400 

NRT 

8280.00 

24.81 

7599 ± 400 

740 ± 103 

8.90 ± 1.30 

J1858-2216 

2.38 

0.39 

820 

GBT 

113.00 

6.63 

qaq +200 

7.64 ± 1.57 

6.80 ± 1.40 

J1902-5105 

1.74 

0.90 

1400 

PKS 

686.00 

22.13 

1 1 79+220 

35.9 ± 2.5 

5.20 ± 0.40 

J1939+2134^ 

1.56 

10.55 

1400 

NRT 

10970.00 

99.54 

3559 ± 350 

139 ±50 

1.30 ±0.50 

J1959+2048^ 

1.61 

0.81 

1400 

NRT 

763.30 

25.49 

2490+1^^ 

126 ± 14 

16.5 ± 1.8 

J2017+0603 

2.90 

0.83 

1400 

NRT 

130.00 

5.93 


98.2 ± 6.3 

75.6 ±4.8 

J2043+171F 

2.38 

0.43 

1400 

AO 

126.50 

6.99 

1700+150 
1 /OU_32o 

100 ±6 

79.2 ± 4.7 

J2047+1053 

4.29 

2.10 

820 

GBT 

105.00 

3.54 

2050!|?o 

31.0 ±6.9 

29.5 ± 6.6 

J205 1-0827^ 

4.51 

1.26 

1380 

WSRT 

54.30 

2.42 

1040!‘“ 

4.37 ± 1.33 

8.10 ±2.50 

J2124-3358^ 

4.93 

1.12 

1400 

NRT 

37.00 

1.83 

300!™ 

3.96 ±0.17 

10.7 ± 0.5 

J2214+3000 

3.12 

1.50 

1400 

NRT 

192.00 

6.63 

1540!™ 

92.9 ± 4.3 

48.4 ± 2.2 

J2215+5135 

2.61 

2.34 

2000 

GBT 

519.00 

12.95 

3010!™ 

128 ± 16 

24.6 ± 3.0 

J2241-5236 

2.19 

0.87 

1400 

NRT 

260.00 

12.24 

Cl O+80 

J1-5_76 

10.5 ± 0.5 

4.00 ± 0.20 

J2302+4442 

5.20 

1.33 

1400 

NRT 

38.20 

1.74 

11VU_23o 

62.2 ± 2.9 

163 ±8 


Notes. 

^ P, E, and 5 lc have been corrected for the Shklovskii effect and Galactic acceleration as in 2PC. 

^ Observatory responsible for generating the radio profile: AO = Arecibo Observatory, GET = Green Bank Telescope, NRT = Nangay Radio Telescope, PKS = 
Parkes Radio Telescope, and WSRT = Westerbork Synthesis Radio Telescope. 


power (E = 4n^IP/P^) values reported in columns 3 and 6 
are corrected for the Shklovskii effect (Shklovskii 1970) and 
Galactic acceleration following 2PC. The observed frequency 
and observatory responsible for generating the radio profile are 
given in columns 4 and 5, respectively. The corrected P val- 
ues are used to calculate the magnetic field strength at the 
light cylinder (5 lc = 47r^((1.5/P)/(c^P^)), for an orthogo- 
nal rotator) given in column 7. The gamma-ray luminosities 


(Ly = ^TtfaGiood^, with Gioo the 2PC gamma-ray energy flux 
from 0.1 to 100 GeV, d the pulsar distance, and /q a beaming 
factor defined in Section 8.1) in column 8 are derived using 
the distance estimates adopted in 2PC. Column 9 gives the effi- 
ciency with which rotational energy is turned into gamma rays 

(r]y = Ly/E). 

Figure 1 presents the LAT-detected MSPs, labeled by model 
class, on a P-P diagram with lines of constant ^lc and E as 
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Table 5 

Simulation and Fitting Parameters for each MSP 


PSR 

P sim 

(ms) 

^sim 

(MHz) 

riy 

(Bins) 

flR 

(Bins) 

by 

(Counts /Bin) 

bR 

(Arb. Units/Bin) 

(Arb. Units) 

Model Class 

J0023+0923 

3.5 

2000 

15 

30 

15.1060 

0.8550 

3.2930 

I 

J0030+0451 

4.5 

1400 

90 

90 

21.5035 

0.1400 

0.6534 

I 

J0034-0534 

1.5 


45 

45 

11.628 

0.680 

0.472 

II 

JOlOl-6422 

2.5 

1400 

30 

60 

9.1766 

1.6400 

13.4850 

I 

J0102+4839 

2.5 

1400 

30 

30 

11.3830 

0.0060 

0.0094 

I 

J021 8+4232 

2.5 

1400 

30 

60 

79.1190 

0.1972 

0.1250 

I 

J0340+4130 

3.5 

1400 

60 

60 

7.5300 

1.6790 

2.8060 

III 

J0437-4715“-'’ 

5.5 

1400 

90 

90 

8.1020 

0.0000 

0.6743 

I 

J0610-2100^ 

3.5 

1400 

20 

20 

8.5540 

1.4460 

2.1335 

I 

J0613-0200^ 

3.5 

1400 

60 

60 

11.5090 

0.2340 

1.6630 

I 

J0614-3329 

3.5 

1400 

90 

180 

28.0340 

3.0710 

1.5770 

I 

J0751+1807 

3.5 

1400 

45 

90 

6.6950 

0.6370 

2.8380 

I 

31024-0719"= 

5.5 

1400 

15 

15 

4.5980 

1.5440 

6.2700 

I 

J1 124-3653 

2.5 

800 

30 

30 

8.0760 

0.0000 

0.0007 

I 

J1 125-5825 

3.5 

1400 

15 

30 

6.3028 

0.0040 

0.0122 

I 

J1231-1411 

3.5 

1400 

90 

90 

28.8400 

0.3990 

0.3540 

I 

J1446-4701 

2.5 

1400 

15 

30 

9.1440 

0.0018 

0.0025 

I 

J1514-4946 

3.5 

1400 

60 

60 

15.3750 

0.4320 

0.6260 

I 

31600-3053^’"= 

3.5 

1400 

30 

30 

6.1387 

0.6525 

6.4656 

I 

31614-2230 

3.5 

1400 

60 

60 

8.9790 

2.1540 

10.0710 

I 

31658-5324"^ 

2.5 

1400 

45 

45 

19.1750 

1.3390 

14.8830 

I 

31713+0747^ 

4.5 

1400 

15 

30 

14.0380 

0.5160 

8.9440 

I 

31741+135F 

3.5 

1400 

15 

30 

2.8650 

0.0660 

0.8280 

III 

31744-1134 

4.5 

1400 

45 

45 

20.2670 

0.6060 

10.4710 

III 

J1747-4036 

1.5 

1400 

15 

15 

7.8620 

3.3660 

64.3190 

I 

31810+1744 

1.5 


45 

45 

18.413 

0.525 

0.708 

II 

J1823-3021A 

1.5 


30 

30 

5.948 

1.924 

1.457 

II 

31858-2216 

2.5 

800 

30 

60 

4.9741 

0.0026 

0.0021 

III 

31902-5105 

1.5 


45 

45 

16.729 

1.878 

21.175 

II 

31939+2134 

1.5 


30 

30 

290.938 

0.000 

0.284 

II 

31959+2048 

1.5 


45 

45 

9.521 

0.691 

0.765 

II 

32017+0603 

2.5 

1400 

60 

180 

9.6560 

0.2380 

0.0926 

I 

32043+1711"= 

2.5 

1400 

60 

60 

12.1462 

0.1135 

0.0505 

I 

32047+1053"^ 

4.5 

800 

20 

60 

11.0960 

0.0031 

0.0096 

I 

32051-0827 

4.5 

1400 

15 

15 

8.1139 

0.4159 

1.0245 

I 

32124-3358 

4.5 

1400 

60 

60 

14.3570 

0.3410 

1.8000 

III 

32214+3000 

3.5 

1400 

45 

90 

20.3950 

2.2670 

1.0987 

III 

32215+5135 

2.5 

2000 

15 

30 

17.9653 

7.5109 

5.0260 

I 

32241-5236 

2.5 

1400 

45 

90 

19.8430 

0.0360 

0.4180 

I 

32302+4442 

5.5 

1400 

90 

90 

10.4260 

1.4130 

2.2020 

I 


Notes. 

^ Polarization suggests presence of a core beam. 

^ Non-standard ctr estimate, one-third the normal value used. 

^ Non-standard ctr estimate, one-half the normal value used. 

^ Non-standard ctr estimate, twice the normal value used. 

^ Non-standard ctr estimate, one-fourth the normal value used. 


well as MSPs from the ATNF pulsar catalog^^ (Manchester et al. 
2005) with no pulsed gamma-ray detection in 2PC. 

Table 5 gives the simulation and fitting parameters for each 
MSP For Class I and II MSPs, we always use a hollow-cone 
radio model and only include a core beam if polarization data 
in the literature suggests there is evidence for such a component 
(noted in the table). For Class II MSPs Vsim is not a parameter of 
the simulations, and thus those entries are empty in Table 5. We 
use events with reconstructed directions within 2° of each MSP 
with spectrally derived probabilities of coming from the pulsar 
^0.2, except for PSRs J1939-I-2134 and J2047-I-1053 for which 
we impose a minimum probability of 0.01 and 0.1, respectively, 
in order to optimize the unweighted counts light curves of these 


http://www.atnf.csiro.au/research/pulsar/psrcat 


two MSPs. From these events, we produced folded light curves 
with Hy fixed- width bins. 

The light curves presented in 2PC use weighted counts (where 
the weight of each photon is the probability that it originated 
from the pulsar based on a spectral and spatial model of the 
region, Kerr 2011) while we use unweighted counts. We made 
this choice as it allows us to use Poisson likelihood for the 
gamma-ray light curves. Poisson likelihood more naturally 
matches sharp peaks, particularly for faint MSPs, than a 
statistic, which is necessary when using binned weighted-counts 
light curves. While 2PC used the weighted H-test value to 
determine the number of bins in the gamma-ray light curves, 
we have adopted a different approach. We only use numbers 
of bins which are an integer divisor of 180 (the number of 
bins in our simulations), and we choose the number for each 
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Table 6 

TPC Model Fit Results for Class I MSPs 


PSR 

-In (C) 

ndof 


a 

(o) 

(0) 

w 

(% 0pc) 

fa 

10023+0923^ 

58.4 

39 

0.567 

38!.'’n 

65^2 

5.0 ± 5.0 


J0030+0451 

581.1 

174 

0.044 

74 ±2 

55!3j 

2.5 ± 2.5 

1 1 Q+O-07 
^•^■^-0.01 

JOlOl-6422^ 

131.7 

84 

0.333 

30+‘2‘ 

84+52 

0.072255 

0 74+0.33 
^•'^-0.07 

10102+4839^ 

102.0 

54 

0.100 

43+35 

^■5-28 

7rj+ll 

'^-16 

0-0-255 

^•^'^-0.01 

J0218+4232 

325.3 

84 

0.033 

25+2 

12 ±2 

0.07*2°5° 

0 1 9+0.01 
^•^^-0.08 

J0437-4715^ 

349.3 

173 

0.022 

35 ± 1 

64 ± 1 

50.0725-50 

'+'•''+'-0.01 

10610-2100^ 

70.8 

34 

0.000 


49+17 

^^-21 


i-ina 

10613-0200^ 

218.8 

113 

0.067 

55+‘j 

43+2 

2.5 ± 2.5 

0.78%% 

J0614-3329 

1312.5 

264 

0.067 

63 ± 1 

84 ± 1 

0.0 ± 2.5 

1.03 ±0.01 

J0751+1807 

206.3 

129 

0.100 

21+2° 

69 ±2 


o -(> o%Z 

11024-0719^ 

50.8 

24 

0.033 

66^65 

70+11 

'^-13 

0.07*2°5° 

0.83%% 

11124-3653^ 

102.9 

54 

0.067 

13+“*/ 

AQ+3 

2 5+2-5 

0 55+2-45 
^•^^-0.13 

11125-5825^ 

40.9 

39 

0.033 

29+f, 

71+10 
' ^-13 


'+'•'^'+'-0.03 

J1231-141F 

835.7 

174 

0.022 

26+3 

69 ± 1 

0.0 ± 2.5 

^•^'+'-0.01 

J1446-470F 

55.6 

39 

0.033 

17+” 

68+5,0 

2 5+^-5 

Q C7+O.6I 
^•^'-0.09 

11514-4946 

212.8 

114 

0.450 

24+’2 

68 ± 1 

0.0 ± 2.5 

0 54+0-04 

^•2’^-0.05 

11600-3053 

84.3 

53 

0.033 

61:‘3° 

37!5„ 


'+'•'^-0.02 

11614-2230 

207.4 

114 

0.167 

00 

o 

00 

0 

1 + 

7.5 di 2.5 

^•'+'^-0.09 

11658-5324^ 

159.0 

84 

0.411 

30+\2 

6975,3 

7 

'•^-7.5 

0-607°o*o8 

11713+0747^ 

44.8 

38 

0.033 

36+5 

68 ±4 

iO-072,-io 

0-687°o5(,^ 

j1747-4036^ 

46.0 

24 

0.167 

iq+71 

35'_18 

807*8°o 

0.07*2°5° 

0.487*o52 

12017+0603^ 

438.5 

234 

0.150 

34+3 

6772j 


0.64%% 

12043+17 IF 

245.8 

114 

0.767 


16% 

7-5752°5 

0-837°oS 

12047+1053^ 

66.0 

74 

0.050 

51+‘; 

717‘i* 

0-072255 

Q c^+0.63 
^•^'-0.05 

12051-0827^ 

48.9 

24 

0.500 

422+47 

^^-22 

AQ+21 

OV_42 

0.07*2°5° 

0-5273o'o*5 

12215+5135 

50.8 

39 

0.900 

18+20 

717*1 

0.07^5 

0.487°oi 

12241-5236^ 

233.0 

129 

0.056 

20 ± 1 

76 ± 1 

2.5 ± 2.5 

0 6T+0-02 

^•^■^-0.03 

12302+4442 

383.4 

174 

0.200 

60+32 

46735 

2.5 ± 2.5 

n 97+0.06 

u.yj_o 11 


Note. ^ Confidence contours in ol and ^ are not simply connected or simple shapes, quoted uncertainties should be used with caution. 



Figure 1. P vs. P, LAT-detected MSPs are labeled by model class as indicated: 
Class I MSPs are black squares. Class II MSPs are light gray (red in the online 
version) upward pointing triangles, and Class III MSPs are dark gray (blue in 
the online version) downward pointing triangles. Gray circles are known radio 
MSPs with no detection in 2PC. Dashed lines show constant Plc while dot- 
dashed lines show constant E, both assuming / = 10^^ g cm^. All but one 
Class II MSP has P <2 ms. 

(A color version of this figure is available in the online journal.) 


MSP, largely by eye, to be as large as possible while preserving 
the significance of structures for the fitting procedure. We have 
also decreased the binning of the radio profiles to be an integer 
divisor of 1 80 bins and to smooth some of the fine structure that 
our geometric models cannot reproduce. 

The gamma-ray background levels were estimated as by = 
{N — S)/ny, where is the total number of events passing the 
selection criteria for each MSP, S' is a proxy for the pulsar signal 
taken to be S = pi with pt the probability that event i is 
associated with the pulsar, and is the number of bins in the 
light curve. The radio background levels (Z?r) were estimated 
from fits of a constant value to the off-peak regions of the radio 
profiles using the estimated radio uncertainties (ctr, derived as 
described in Section 5 with exceptions noted in Table 5) divided 
by the number of radio bins (^r). 

For the light curve fits, we matched the observed period to 
the closest simulated period (Psim) except for Class II MSPs 
for which we always use Psim = 1.5 ms. For all but one of the 
Class II MSPs, this choice of Psim matches the observed period 
well, the exception being PSR J1 823— 3021 A with a period 
of 5.44 ms. However, the period enters the simulations most 
strongly through the size of the polar cap, which means that 
using a shorter period will, at most, overestimate the size of any 
predicted off-pulse region. For Class I and III MSPs, we also 
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Table 7 

OG Model Fit Results for Class I MSPs 


PSR 

-In (£) 

ndof 


a 

(o) 

(o) 

^acc 

(% 0Pc) 

^em 

(% 0Pc) 

fa 

J0023+0923^ 

60.7 

38 

0.567 

67 ±7 

24+10 

0.0/;/ 

0.0 ± 2.5 

0-81±53 

J0030+0451 

506.9 

173 

0.011 

00 

00 

68 ± 1 

2.5 ± 2.5 

0.0 ± 2.5 

0-98T03 

JOlOl-6422^ 

135.7 

83 

0.833 

90 ±21 

35/; 

5.0 ± 5.0 

2.5 ± 2.5 

0 QQ+0.08 

J0102+4839^ 

97.8 

53 

0.133 

60+30 

76/ 

0.0/« 

0.0 ± 2.5 

0 73+0-39 

J021 8+4232^ 

2355.2 

83 

0.167 

45+33 

^^-35 

67+23 

0/_47 

0-0-2" 

0.0 ± 2.5 

0.52/“ 

10437-4715^ 

340.1 

172 

0.022 

76 ± 1 

46 ± 1 

10.0 ± 2.5 

2.5 ± 2.5 

0.82// 

J0610-2100^ 

79.6 

33 

0.500 

63+13 

02)_4o 

89/21 


0.0 ± 2.5 

1.UZ_o74 

J0613-0200^ 

226.6 

112 

0.067 

60 ± 1 

45 ± 1 

5.0 ± 2.5 

2.5 ± 2.5 

0 

'^•'^'-0.34 

J0614-3329 

1487.9 

263 

0.089 

00 
^ + 

88 ± 1 

0.0 ± 2.5 

0.0 ± 2.5 

0.98 ± 0.01 

J0751+1807^ 

207.4 

128 

0.033 

59+^2 

72/i 

0-0-2" 

0.0 ± 2.5 

0 55+0.59 

u.jj_0 22 

J1024-0719 

48.2 

23 

0.033 

66!.?8 

70 + 11 
'■^-37 


O.O-40 

0 64+0-60 

J1 124-3653^ 

108.5 

53 

0.067 


73/b 

5.0 ± 5.0 

2.5 ± 2.5 

A 04+1.66 

J1 125-5825^ 

40.2 

38 

0.033 

00 

6O/38 

2 5+^-5 

0.0 ± 2.5 

1 1 7+0.02 
'-0.88 

J1231-141F 

696.8 

173 

0.056 

88 ± 1 

67 ± 1 

2.5 ± 2.5 

0.0 ± 2.5 

0 QQ+0.12 

J1446-470F 

51.7 

38 

0.033 

oo 

o 

00 

25/; 

5.0 ± 5.0 

0.0 ± 2.5 

^■'^- 0.36 

J1514-4946^ 

221.6 

113 

0.450 

25^7 

68 ± 1 

0.0/;/ 

0.0 ± 2.5 

A oc+0.50 
'^•■^^-0.06 

J1600-3053 

84.9 

52 

0.033 

65+3 

28/j 

O.O/o 

0.0 ± 2.5 

0.76/38 

J1614-2230 

204.1 

113 

0.150 

o^ 

O 

fsl “Ci 
+ 1 
OO 
OO 

2-5/“5 

0.0 ± 2.5 

0.88/“ 

J1658-5324 

165.6 

83 

0.411 

78+3, 

23/5 


2.5 ± 2.5 

0.84// 

J1713+0747^ 

50.3 

37 

0.033 

36+'/ 

65/“' 

0.0/;/ 

0.0 ± 2.5 

A 40+0.71 

11747-4036^ 

45.3 

23 

0.167 

10+71 

^^-13 


0.0/;/ 

0.0 ± 2.5 

0 32+1-H^ 

J2017+0603^ 

409.3 

233 

0.150 

23 ± 1 

74 ± 1 

5.0 ± 2.5 

2.5 ± 2.5 

0 30+0-16 

J2043+171F 

303.6 

113 

0.733 

CO +8 
^■^-26 

79/4 

O.O/o 

0.0 ± 2.5 

0.90/3/g 

J2047+1053^ 

67.6 

73 

0.050 

51// 

71/; 

0.0 ± 2.5 

0.0 ± 2.5 

O.38//5 

J205 1-0827 

44.4 

23 

0.500 

68+33 

42+48 

^^-14 

0.0/;/ 

0.0 ± 2.5 

1-08// 

J2215+5135 

48.8 

38 

0.900 

18// 

71/; 


0.0 ± 2.5 

A aa+0.41 
'^•^^- 0.01 

J2241-5236 

287.7 

128 

0.056 

19 ± 1 

75 ± 1 

0.0 ± 2.5 

0.0 ± 2.5 

0.44 ± 0.01 

J2302+4442^ 

390.2 

173 

0.189 

65/ 

37/ 

0.0/« 

0.0 ± 2.5 

A qo+0.16 
'^•^°-0.03 


Note. ^ Confidence contours in ol and ^ are not simply connected or simple shapes, quoted uncertainties should be used with caution. 


Table 8 

alTPC Model Fit Results for Class II MSPs 


PSR 

-In (C) 

ndof 

i>/+ 

a 

( 0 ) 

( 0 ) 

Wy 

(% 0 pc) 

Wr 

(% 0 pc) 

+'^max 

(^LC) 

dR 

'^min 

(^LC) 

dR 

^max 

(^LC) 

fa 

J0034-0534^ 

J1810+1744^ 

J1823-3021A^ 

J1902-5105^ 

J1939+2134^ 

J1959+2048^ 

130.9 

150.0 

83.4 

164.2 

137.3 

130.4 

80 

80 

50 

80 

50 

80 

0.567 

0.611 

0.633 

0.189 

0.433 

0.522 

23/3 
82/, 0 

46+20 

^°-38 

10/ 

88/ 

56+“ 

2’n-i6 

69 ±3 

16+^ 

11^-10 

52+17 

•^^-28 

73/ 

88/ 

85/ 

2 5+^-^ 
^•2*-2.5 

5.0 ±5.0 

2 5+^-5 

0 - 0 - 2 / 

2 5+^-^ 
^•2*-2.5 

2 5+^-5 
0.0+^2/0 
0.0+^2/0 
0.0 ± 2.5 
0.0+^2/0 

A ac+0.25 

0.70/“ 
0.85/35 
0.95 ± 0.25 
0.85/35 

0-95/1 

0.85/05 

0.14/1 

0.65/30 

0-80/'0 

A 1 4 + 0.81 

1 oo+n-^n 

l-'^'^-0.30 
0.70 ± 0.50 
0.70 ± 0.50 

A ac+0.05 

0.95 ±0.10 

A ac+0.25 

A 4^+0.59 

1-^1 -0.56 
A 04 + 0.95 
l^-^^-0.36 

A 49+1.71 

l-l^^-0.03 
0 RQ+1-25 


Note. ^ Confidence contours in oi and ^ are not simply connected or simple shapes, quoted uncertainties should be used with caution. 


matched the simulated radio frequency (Vsim) to the observed 
frequency as closely as possible since this affects the assumed 
emission height and cone/core size. For Class II MSPs, we do 
not specify Vsim as we assume all frequencies come from an 
extended region in the magnetosphere. 

7. RESULTS 

The best-fit parameters for Class I MSPs are given in Tables 6 
and 7; for Class II MSPs in Tables 8, 9, and 10; and for Class III 
MSPs in Table 11. In each table, we give the phase of the 


magnetic axis (O^) predicted from the models. The observed 
and best-fit light curves for all 2PC MSPs are provided in 
Appendix A along with information about the initial radio 
and gamma-ray discoveries and a comparison of the best-fit 
geometries to other estimates in the literature, when they exist, 
for each MSP. 

Sample confidence contours are shown in Appendix B. 
All confidence contours are included in the auxiliary online 
material. For contours that are not simply connected, the 
uncertainties on a and ^ are estimated only from the region 
that contains the best-fit geometry. As such, uncertainties for 
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Table 9 

alOG Model Fit Results for Class II MSPs 


PSR 

- ln(/:) ndof 


a 


^y,a.cc 

n’y,em n>R,acc 

n'R,em 

TfV 

+'^max 

nR 

+'^min 

nR 

+'^max 

/q 




(o) 

(o) 

(% 0Pc) 

(% 0pc) (% 0pc) (% 0pc) 

(/?LC) 

(^LC) 

(^LC) 


10034-0534^ 

137.9 78 

0.589 

22+^^ 

zz_i7 

75+‘| 

5.0 ±5.0 

2.5 ±2.5 5.0 ±5.0 2.5 ± 2.5 

0 95+0-25 

0 QO+0-05 

^•^'9_o.76 

1 05+0-^5 

0 24+0-^! 
^•^^-0.06 

11810+1744'“ 

147.8 78 

0.167 

o 

00 

2U'i4 

2 

^•^-2.5 

0.0 ±2.5 5.0 ±5.0 2.5 ±2.5 

0.95 ± 0.25 0.25+'Jj’|®i 

1 1Q+0.10 

0 90+0-15 

'5-5'+'-0.72 

J1823-3021A 

81.2 48 

0.467 

42!^36 

28^*65 

5.0 ±5.0 

2.5 ±2.5 5.0 ±5.0 2.5 ± 2.5 

1 1Q+0.10 

0.95+0j|, 

1 00+0-20 
^•^^-0.45 

0 45+2.23 
^•^^-0.27 

11902-5105^ 

159.2 78 

0.189 

10 ±2 

73 ± 1 

0-0-2.5 

0.0 ± 2.5 O.O^z^s 

0.0 ± 2.5 

0.95+«i‘0 

0-80+«o‘a^6 

0 95+0.10 

0 30+0-22 
'5-5'5-0.10 

11939+2134^ 

135.7 48 

0.433 

72 ± 1 

85 ± 1 

0.0+'2®5° 

0.0 ± 2.5 2.5+” 

0.0 ± 2.5 

0.80!°,% 

0 90+0-05 

1.05 ±0.15 

1-01±63 

J1959+2048 

131.8 78 

0.544 

OZ_42 

87-13 


0.0 ± 2.5 0.Q%§ 

' 0-0-2" 

1.00%% 

0 1 4+0.81 

1 OO+0-20 
1-'9'9-0.40 

0-80%% 

Note. ^ Confidence contours in a 

and ^ are 

not simply connected or simple shapes, quoted uncertainties should be used with caution. 
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laSG Model Fit Results for Class II MSPs 
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-In (C) 
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a 


Wy 

Wr 

^in 

^out 

Sf 

/q 





(o) 

(o) 

(% 0Pc) 

(% 0pc) 

(^NS) 

(^NS) 

(^NS) 


10034-0534^ 

133.3 

80 

0.033 

22+9 

74+4 f ) 0 + 5.0 

/^_34 ^-^- 2.5 

0.0+5, 0 


0.9 ± 0.6 

1 O + 0.2 
^•^- 0.6 

0.86+Oj« 

J1810+1744'' 

150.7 

80 

0.033 

c+81 

^-3 

29+47 5 0 + 2.5 

Z ^-27 ^-^- 5.0 

2.5 ± 2.5 


0 9 +^-^ 

1 q + 0.2 
^•^- 0.6 


J1823-3021A^ 

87.7 

50 

0.567 

78;!:‘|' 

45+6 5 

^ 0_22 ^•' 4 - 5.0 

0.0+5, 0 

0 1 + 0-0 

0 3 +^-^ 
- 5 - 0.3 

1 ^+ 0.6 
5 -^- 0.2 

0 60 +^ 53 
U . DU _ o 28 

11902-5105^ 

166.1 

80 

0.167 

56 ± 8 16+*, 

0 

U . U _25 

5-0-¥o 

0 3 + 0-5 
- 5 - 0.2 

2 1 + 0-5 

1.6 ±0.4 

0 92 + 0-62 
'+•^^- 0.39 

11939+2134^ 

154.3 

50 

0.833 

33 ± 12 41+1 

, U . U _25 

2.5 ± 2.5 

0 2 + 0-5 

0 3 + 0-6 
- 5 - 0.3 

1 4 + 0.6 
5 -^- 0.2 

1 Cl + 0.40 
^•^^- 0.55 

11959+2048 

125.5 

80 

0.433 

46^38 

35% 0.0+50 

2.5 ± 2.5 

0 3 + 0-2 

u .:^_02 

2 1 + 0-5 

1.6 ±0.4 

1 63+^-^^ 

i. 0.5-1 17 


Note. ^ Confidence contours in ol and ^ are not simply connected or simple shapes, quoted uncertainties should be used with caution. 


Table 11 

PSPC Model Fit Results for Class III MSPs 


PSR 

-In (£) 

ndof 


a 


/q 





(o) 

(o) 


10340+4130^ 

197.4 

115 

0.050 

43+'“ 

73 ±7 

1 ^71^+0.20 
^•^2)_o29 

J1741+135F 

33.4 

40 

0.300 

46 %, 

00 

o 

1 06+^-^^ 
^•'50-0.07 

11744-1134^ 

174.8 

85 

0.033 

51% 

85+3 

1 04+ 1-00 
1-5^^-0.53 

11858-2216 

93.3 

85 

0.067 

42 % 

74+* 

1 i^+0.23 

12124-3358^ 

286.2 

115 

0.983 

19 ± 1 

25+‘*2 

0 46+0.50 

H-^11-0.02 

12214+3000 

207.2 

130 

0.056 

59+‘3 

86 ±2 

1.05 ±0.02 


Note. ^ Confidence contours in ol and ^ are not simply connected or simple 
shapes, quoted uncertainties should be used with caution. 


fits that are flagged as not simply connected should be used with 
care and, when possible, the confidence contours should be 
used directly. For some of the brightest MSPs in our sample, the 
resulting confidence contours, even after rescaling as discussed 
in Section 5, suggest uncertainties on the order of a few 
degrees. For these MSPs, we remind readers of the systematic 
uncertainties discussed in Section 5.1. 

The light curve fitting is only done for o' c [1°, 90° ] and 
^ G [0°, 90°). When using geometric models, the light curves 
produced beyond a — 90° or ^ = 90° are the same but shifted 
by 0.5 in phase, and those with a and ^ both >90° are exactly 
the same. This arises due to the assumption that both poles have 
the same emission properties. Thus, while the uncertainties we 
quote are constrained to the region with a and ^ both ^90°, 
there are identical confidence regions when reflected across the 
90° boundaries. 

Sample skymaps of simulated gamma-ray and radio emission 
(constructed as detailed in Section 4), corresponding to the best- 
fit models of the given MSP, are shown in Appendix C. Emission 


maps for the best-fit models of all MSPs are included in the 
auxiliary online material. 

The gamma-ray and radio light curves of several 2PC MSPs 
have been modeled previously using similar models and meth- 
ods. Venter et al. (2009) modeled the light curves of the first 
eight MSPs detected with the LAT, choosing the best geome- 
tries by eye with steps of 5° in both a and This approach 
put more weight on matching the gamma-ray light curves. Our 
results agree well for PSRs J0030-I-045 1 , J0437— 4715 (the TPC 
fit), and J1744— 1 1 34. For the other MSPs they fit, our likelihood 
method either prefers a different maximum (see the discussion 
of confidence contours with multiple maxima in Appendix B) 
or the best fit is substantially different due to increased statistics 
and/or more emphasis on matching the radio light curve. 

The light curves of PSR J0034— 0534 have been modeled 
by Abdo et al. (2010b) using by-eye fits with alOG and alTPC 
models and by Venter et al. (2012). The latter authors used a 
similar likelihood fitting technique for alOG and alTPC models 
and by-eye fits for the laSG model. In both cases, the alOG and 
alTPC results agree with the best-fit parameters presented in 
Tables 8 and 9. However, our laSG best-fit geometry does not 
agree with that reported by Venter et al. (2012), likely because 
their by-eye fits did not put emphasis on the low-level off-peak 
emission. 

Both previous studies used the 320 MHz radio profile and 
not the 1400 MHz profile we use here. At lower frequency, the 
peaks are at roughly equal heights whereas at 1400 MHz the 
peak near phase ~0.8 is significantly lower than the other peak. 
Compared to the best-fit parameters of Venter et al. (2012), we 
do And slightly higher values of and though the 

differences are not significant when considering the estimated 
uncertainties. This is opposite of what one would expect from 
radius-to-frequency mapping studies (e.g., Cordes 1978) and 
may point to additional differences between the radio emission 
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observed in Class II MSPs and the single-altitude cone and/or 
core assumption. 

Cognard et al. (2011) fit the light curves of PSRs 12017-^0603 
and 12302-^4442 with OG and TPC gamma-ray models, a 
hollow-cone beam, and a Markov chain Monte Carlo maximum 
likelihood procedure. Our results agree well with theirs for 
PSR J2302-I-4442, but we find somewhat larger a values for 
PSR 12017-^0603. 

Light curve fits for PSR J1823— 3021A were reported by 
Freire et al. (2011), in the supplementary online material, using 
the alTPC and alOG models and a similar maximum likelihood 
technique. These fits were done above 500 MeV due to issues 
with the background level estimation. Their best-fit values agree 
with ours within the quoted uncertainties, though our values of 
are noticeably greater than theirs, for both models. 

Guillemot et al. (2012b) reported fits to the light curves of 
PSRs 11939-^2134 (B 1937-^21) and 11959-^2048 (B 1957-^20) 
using alOG and alTPC models and a similar maximum like- 
lihood procedure. Venter et al. (2012) provided more details 
on these fits as well as by-eye fits using the laSG model. For 
PSR J1939-I-2134, the fits agree within uncertainties, but our 
alTPC best-fit geometry agrees better with those from fitting ra- 
dio polarization data (also presented by Guillemot et al. 2012b). 
For PSR 11959-^2048, our alOG and alTPC fits agree with theirs, 
within the estimated uncertainties, while our best-fit laSG model 
has a slightly different geometry (though still within uncertain- 
ties). Their fits used the 350 MHz radio profile and not the 
1400 MHz profile. The lower frequency profile only has two 
peaks, both of which are matched to a peak in the gamma-ray 
light curve. The 1400 MHz profile has an additional peak that 
does not have a corresponding feature in the gamma-ray light 
curve. This additional peak was not fit by our models, but its 
position in phase may indicate that it is from either a core or 
conal component. 

Kerr et al. (2012) modeled the gamma-ray and radio light 
curves of PSR JOlOl— 6422 using TPC and OG gamma-ray 
models with a hollow-cone radio beam. Our fits are compatible 
with theirs, but the models have difficulties producing both the 
gamma-ray and radio light curves. This MSP may require a 
more complex radio model. 

The gamma-ray and radio light curves of PSR 12043-^1711 
were fit by Guillemot et al. (2012a) using TPC and OG gamma- 
ray models with a hollow-cone radio beam and the same 
maximum likelihood procedure. Our best-fit results for this MSP 
are in agreement with theirs. 

Du et al. (2010) modeled only the gamma-ray light curves of 
PSRs 10030-^0451, 10218-^4232, and J0437-4715 using data 
from the first LAT pulsar catalog (Abdo et al. 2010d) and an 
annular gap model, assuming that the magnetic axes of MSPs 
should have moved toward alignment with the spin axis. For 
PSR 10030-^0451, our fits do not agree well with theirs, finding 
larger a values consistent with X-ray light curve modeling, and 
their model overpredicts the observed level of off-peak emission. 
For PSR J021 8-1-4232, our OG fit is consistent with theirs; how- 
ever, our TPC best-fit geometry has a substantially lower ^ value, 
necessary to match the radio profile. For PSR J0437— 4715, our 
TPC fit finds a similar a with somewhat larger ^ . Our OG fit 
finds a significantly larger a and slightly lower ^ . 

Du et al. (2013) modeled the gamma-ray and radio light 
curves of PSRs J0034-0534, JOlOl-6422, and J0437-4715 
using an annular gap model. Their fits were done by eye but 
guided by the assumption that the inclination angles should 
be small (consistent with the possible trend of decreasing a 
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with age for non-recycled pulsars, e.g.. Young et al. 2010) and 
by geometric estimates from other wavelengths when possible. 
Their geometry for PSR J0034— 0534 is in agreement with 
all of our fits, though their model does not produce any off- 
peak emission, which is interesting given that the annular 
gap emission geometry is different from our models. For 
PSR JO 10 1—6422 their models reproduce the data well and 
neither of our best-fit geometries agree with theirs. Finally, 
our best-fit results do not agree with their quoted geometry for 
PSR J0437— 4715 (note that they fixed ^ a priori to agree with 
the estimated orbital inclination angle, van Straten et al. 2001; 
Hotan et al. 2006). Bogdanov et al. (2007) found it necessary 
to use a magnetic dipole that was offset from the center of the 
PC to fit the X-ray light curve of this MSP, a possibility we 
plan to explore in the future. In comparing our fits to those of 
both Du et al. (2010) and Du et al. (2013), it is important to 
note that their simulations used a static dipole magnetic field 
configuration. It is unclear exactly how including the effects of 
rotation will change their best-fit geometries. 

Petri (201 1) modeled the observed gamma-ray and radio light 
curves of eight pulsars (three of which were MSPs) using a 
striped wind model for the HE emission and a polar cap beam 
for the radio profile. For PSR 10030-^045 1 our fits agree well 
with those of Petri (2011) in that both a and ^ should be >50°, 
but we find o' > ^, the opposite of what he finds. Additionally, 
we note that he had to increase the size of the polar cap emitting 
region to a significant portion of the entire neutron star surface 
to match the radio peaks and his predicted radio to gamma-ray 
phase lag is too large by 0.1 in phase. For PSR 10218-^4132, 
our OG fit results agree well with his geometry, though his 
model radio peaks are not broad enough and the double-peaked 
structure of the main component is not reproduced. Our best-fit 
results for PSR J0437— 4715 do not agree well with his model 
geometry (except for the ^ value of our OG fit) and while he 
does reproduce the basic light curve structures well the predicted 
phase lag is too large. 

Espinoza et al. (2013) identified PSRs 10218-^4232, 
J2017-F0603, J2124-3358, and 12302-^4442 as W-type, sug- 
gesting that these were aligned rotators with a close to 0°. For 
PSRs J0218-F4232, 12017-^0603, and J2124-3358 the best-fit 
geometries have values of a between 20° and 30° (excluding 
the OG fit for PSR 10218-^4232 which does not reproduce the 
radio profile well), supporting the aligned rotator hypothesis. 
For PSR J2302-I-4232, both the TPC and OG models find a 
near 60°. While neither exactly reproduces the radio profile (see 
Appendix A) the gamma-ray light curve is reproduced well and 
Johnson (201 1) showed that, for the same a and ^ , increasing the 
altitude of the simulated radio emission leads to a better fit of the 
radio data. PSR J1 124— 3653 would be identified as a W-type 
MSP and, for both the TPC and OG model fits, is found with a 
low O', between 10° and 20°, which may support the separation of 
W-type and N-type MSPs based on the characteristics of their 
radio profiles. However, there are several Class I MSPs with 
low values of o', in some fits, that would likely be, or have been, 
classified as N-type, namely, PSRs J0751-hl807, J1446— 4701, 
J1514-4946, J2215-F5135, and J2241-5236. If our best-fit a 
values are reliable, this brings into question the morphological 
distinction between the N-type and W-type MSPs of Espinoza 
etal. (2013). 

7.7. Best-fit Parameter Distributions 

Eigure 2 presents the — ln(£) differences between the best 
and next-best models for all Class I and II MSPs. We say that 
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J2302+4442 
J2241-5236 
J2215+5135 
J205 1-0827 
J2047+1053 
J2043+1711 
J2017+0603 
J1959+2048 
J1939+2134 
J1902-5105 
J1823-3021A 
J1810+1744 
J1747-4036 
J1713+0747 
J1658-5324 
J1614-2230 
J1600-3053 
J15 14-4946 
J1446-4701 
J1231-1411 
J1 125-5825 
J1 124-3653 
J1024-0719 
J0751+1807 
J06 14-3329 
J0613-0200 
J0610-2100 
J0437-4715 
J021 8+4232 
J0102+4839 
JOlOl-6422 
J0034-0534 
J0030+0451 
J0023+0923 



10 '^ 1 10 10 ^ 10 ^ 

A -ln(likelihood) 

Figure 2. Difference between the smallest — ln(/^) value and the next best for all Class I and II MSPs. The symbol for each MSP denotes which fit is best as indicated 
in the legend: TPC, light gray (pink in the online version) filled squares; OG, dark gray (green in the online version) hlled upward pointing triangles; alTPC, light gray 
(pink in the online version) open squares; alOG, dark gray (green in the online version) open upward pointing triangles; and laSG, black (red in the online version) 
open circles. The vertical dashed line indicates a — ln(C) difference of 15, our threshold for a significant preference of a particular model over the others. 

(A color version of this hgure is available in the online journal.) 



CO 


Figure 3. Best-fit (a, pairs for all 40 2PC MSPs. The markers indicate which 
model is best as follows: TPC, light gray (pink in the online version) hlled 
squares; OG, dark gray (green in the online version) hlled upward pointing 
triangles; PSPC, black (blue in the online version) hlled downward pointing 
triangles; alTPC, light gray (pink in the online version) open squares; alOG, 
dark gray (green in the online version) open upward pointing triangles; and 
laSG, black (red in the online version) open circles. Only the geometry for the 
ht with the smallest — ln(£) is shown for each MSP. The dashed diagonal line 
indicates the magnetic axis with a = ^. 

(A color version of this hgure is available in the online journal.) 

one model is significantly favored over the other for only seven 
MSPs; those with A(— ln(£)) ^15 are indicated by the dashed 
vertical line. 

Figure 3 presents the best-fit (a, pairs for each MSP, with 
only the values corresponding to the fit with the smallest — ln(£) 


shown for each MSP, even if the log-likelihood difference is 
<15. There is a clear preference for ^ near 90°. This is expected 
if pulsar spin axes are distributed randomly with respect to the 
Earth line of sight and for outer-magnetospheric gamma-ray 
models where the brightest emission is predicted near the spin 
equator. However, the best-fit a values seem to be uniformly 
distributed between 0° and 90°, which may have implications 
for MSP evolution after the recycling phase (see Section 8.2). 

This contrasts with best-fit a > 40° values found for non- 
recycled gamma-ray pulsars by Pierbattista et al. (2014). While 
the a values are spread out in the best-fit TPC models, the best- 
fit OG models generally seem to require a > 60° . This reflects 
the fact that emission in the OG models fills a relatively smaller 
part of (ckf, O phase space, so that pulsars are only visible over a 
restricted range of o'. In the case of the PSPC models, most of the 
best-fit a values lie in the range from 40° to 60°. This gives the 
optimal level of off-peak emission and radio peak multiplicity. 

Figure 4 presents the best-fit accelerating gap widths for 
Class I and II MSPs; only the values corresponding to the fit with 
the smallest — ln(£) are shown. Values from TPC and alTPC 
fits are grouped in the same histogram, as are values from OG 
and alOG fits. The uncertainties in the tables suggest that, in 
general, we cannot constrain the gap width parameters well, but 
it is clear that the light curve shapes prefer smaller widths. Gaps 
with width 0 (i.e., emission only along one field line for each 
/ovc) are unphysical and suggest that the actual best-fit value is 
smaller than our resolution of 2.5% 0pc. 

Figure 5 plots the best-fit emission altitudes for alTPC and 
alOG fits of all Class II MSPs, regardless of which fit is preferred 
by the likelihood in order to better demonstrate the behavior for 
both models. The minimum radio emission altitudes for alOG 
fits in Figure 5 are plotted as either the best-fit value or the 
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Figure 4. Best-fit accelerating gap widths for Class I and II MSPs. Only the 
value for the fit with the smallest — ln(£) is shown for each MSP. TPC and 
alTPC fit values are shown in the light gray filled histogram (pink in the online 
version, filled with downward diagonal lines from left to right). OG and alOG 
fit values are shown in the dark gray (green in the online version) histogram 
hlled with upward diagonal lines from left to right. Values from laSG hts are 
shown in the black (red in the online version) histogram hlled with horizontal 
lines. 

(A color version of this hgure is available in the online journal.) 

average altitude of the NCS for the best-fit a, whichever is 
larger. 

The radio and gamma-ray emission regions typically have 
significant overlap, with the radio emission often concentrated 
toward the upper edge of the gamma-ray emission region. These 
altitudes are significantly greater than typically assumed for 
other radio emission models (e.g., the cone and core beams we 
use for Class I and III MSPs; see Equation (2)). With the current 
statistics and use of the 95% confidence level uncertainties we 
cannot constrain these altitudes well. 

8. DISCUSSION AND CONCLUSIONS 

The gamma-ray MSPs in 2PC (representing approximately 
one-third of the known gamma-ray pulsar population) have 
a greater variety of light curve morphology than the young 
gamma-ray pulsars. While most MSP light curves are similar 
to those of young pulsars, showing the correlation of increasing 
radio lag with decreasing gamma-ray peak separation (Abdo 
et al. 2010d, 2013), a much larger number show phase alignment 
of radio and gamma-ray peaks and a new phenomenon, gamma- 
ray peaks leading the radio peaks, not observed in the young 
pulsar population. We find that the MSP light curves resembling 
those of young gamma-ray pulsars are indeed well fit by the 
narrow gap models that best describe young pulsars (Romani & 
Watters 2010; Pierbattista et al. 2014) and by standard radio core 
and cone geometry (Rankin 1983; Arzoumanian et al. 2002). 

We have successfully fit the light curves showing aligned 
gamma-ray and radio peaks with the same models but assuming 
that the radio peaks are also caustics, and fitting for minimum 
and maximum radius of emission. This implies that the radio 
emission may occur at high enough altitude and may be extended 
enough along the last open field lines, that the emission is 
also strongly distorted into caustic patterns. Manchester (2005) 
and Ravi et al. (2010) have suggested that similar models are 



J0034-0534 J1810+1744 J1823-3021A J1902-5105 J1939+2134 J1959+2048 

Figure 5. Best-fit radio and gamma-ray emission altitudes from alTPC and alOG 
fits of all Class II MSPs. The filled areas show the extent of the emission regions. 
The asterisk points show the best-ht values and corresponding uncertainties. The 
horizontal dashed line indicates the neutron star surface at 0. 14 /?lc for a 1 .5 ms 
spin period. The minimum altitudes of the gamma-ray emission regions for 
alOG models are set to the approximate minimum altitude of the NCS for the 
best-fit O'. 

(A color version of this hgure is available in the online journal.) 

needed for the radio emission of all gamma-ray pulsars. One 
prediction of this model, that the caustic peaks should have low 
linear polarization, seems to be borne out by radio observations 
of these pulsars showing that most in fact have little or no 
polarization (also noted by Espinoza et al. 2013). We have fit 
the MSP light curves in which the gamma-ray peaks lead the 
radio peak(s) with PSPC models where gamma-ray emission 
occurs over the whole open-field volume and we are viewing 
the high-altitude portion that is highly aberrated to earlier phase. 

The best fits provide estimates of a and as well as emission 
gap width, for each model. Additionally, the best fits for the 
altitude-limited models provide minimum and maximum radii 
for the radio and maximum radius for the gamma-ray emission. 
We find best-fit a values are nearly evenly distributed over a 
large range between 10° and 90°. On the other hand, our best-fit 
^ values are concentrated mostly at high values between 70° 
and 90°, due to the caustics being brighter at large ^ and the 
greater probability of viewing at large angles. 

The large a range partly reflects the wide radio beam size and 
smaller magnetospheres of MSPs, allowing both the radio and 
gamma-ray beams to be visible over a larger region of phase 
space than for longer-period pulsars. Thus, a larger range of 
impact parameter ^ is allowed, restricting the number of radio- 
quiet MSPs to very few. Indeed, there are no radio-quiet MSPs 
in 2PC, while half of the young gamma-ray pulsar population 
is radio-quiet. This larger range of phase space accessible to 
observations may also partly explain the increased diversity in 
light curve profiles we observe for MSPs. On the other hand, 
the large a range for MSPs could be a result of their spin- 
up evolution (see Section 8.2). Their relatively weak magnetic 
fields compared to young gamma-ray pulsars require some 
suppression, possibly by flux burying (Romani 1990) or crustal 
plate tectonics (Ruderman 1991), during their accretion and 
spin-up phase. Such processes may be accompanied by changes 
in a. 

We find that the classic radio core/cone geometry, with 
aberration and retardation included, is able to reasonably fit 
the properties of most MSP light curves, such as peak widths. 
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Figure 6. Histogram of /q values, regardless of which ht is preferred by the 
likelihood, for the TPC and alTPC models (light gray filled histogram with 
downward diagonal lines from left to right, pink in the online version), OG and 
alOG models (dark gray histogram with upward diagonal lines from left to right, 
green in the online version), PSPC model (histogram filled with vertical lines, 
blue in the online version), and laSG model (histogram filled with horizontal 
lines, red in the online version). 

(A color version of this figure is available in the online journal.) 


number, and gamma/radio lag. However, there is indication 
for some sources that the radio emission altitude we have 
assumed is too low, possibly suggesting that standard radius- 
to-frequency mapping is not applicable to MSPs (suggested by 
Kramer et al. 1999). In some cases where the radio light curves 
have multiple peaks, such as PSRs J0613— 0200, J1658— 5324, 
and J2302+4442; both the width and radio lags are too small 
in our models. It may be interesting to explore models with 
higher-altitude radio cones and with partially filled cones (Lyne 
& Manchester 1988), which Kramer et al. (1999) suggested 
could explain the apparently abnormal profile development with 
frequency of radio MSPs. Additionally, turning on different sets 
of field lines in a partially filled cone will lead to different 
predicted phase lags. 

Our finding that the majority of MSP light curves are best fit 
by the TPC or OG models with narrow emission gaps was not 
expected prior to the launch of Fermi. The formation of narrow 
accelerator gaps in these models requires sufficient multiplicity 
of electron-positron pairs to screen the parallel electric field 
throughout the rest of the open-field region outside the gaps. 
The standard versions of these models assume dipole magnetic 
fields and the death lines for pair cascades lie above the P of 
most MSPs (Harding et al. 2005; Zhang et al. 2004), suggesting 
they should be pair- starved. We find that only a small fraction 
of MSP light curves are best fit with PSPC models, indicating 
that the pair death lines predicted by the standard models are 
too high. One idea for lowering the pair death lines of MSPs 
is to invoke non-dipolar fields near the surface (Harding & 
Muslimov 2011a, 2011b; Takata et al. 2012) that can produce 
off-set dipoles or polar caps, and larger parallel electric fields 
and radii of curvature. Additional pair-production mechanisms 
than were traditionally assumed have also been explored, such as 
magnetic pair production at low altitudes in OG models (Takata 
et al. 2010) and photon-photon pair production near the surface 
in polar cap models (Zhang & Qiao 1998; Harding et al. 2002). 



Figure 7. Efficiencies corrected for best- fit /q values, excluding the high 
efficiency outlier PSR J061 0—2100 for visualization purposes. Only the value 
for the fit with the smallest — ln(£) is shown for each MSP. TPC and alTPC fit 
values are shown in the light gray filled histogram (pink in the online version) 
filled with downward diagonal lines from left to right. OG and alOG fit values 
are shown in the dark gray (green in the online version) histogram filled with 
upward diagonal lines from left to right. Values from laSG fits are shown in the 
black (red in the online version) histogram filled with horizontal lines from the 
PSPC fits in the black (blue in the online version) histogram filled with vertical 
lines. 


(A color version of this figure is available in the online journal.) 


8.1. Efficiencies 

In Tables 6 through 1 1 , for each MSP and model, are estimates 
of the beaming correction factor /q (Equation (4) of Watters 
et al. 2009, reproduced in Equation (5)), which accounts for 
the fact that MSP emission is not isotropic. The /q values for 
each fit are shown in Eigure 6; here, we do not show only the 
best fits in order to better demonstrate the range of /q values 
corresponding to realistic light curve shapes for each model. We 
estimate /q by summing the simulated emission (Fsim) for all 
viewing angles at the best-fit a and dividing by twice the sum 
of the simulated emission for the best-fit a and ^ . 


/q(«, O = 


/ / fsim(«, (p) sin(C)dCd(p 
2 f Fsim(a, ^ 4>)d(f) 


(5) 


Uncertainties for /q are estimated from the range of values 
found when calculating /q at each {a, i;) pair within the 95% 
confidence-level contours. The /q values for gamma-ray pulsars 
are typically close to and/or consistent with 1 (e.g., Watters et al. 
2009; Venter et al. 2009; Watters & Romani 2011) motivating 
the use of this value in 2PC when reporting Ly . There are cases 
where this assumption leads to efficiencies greater than 100% 
(see column 10 of Table 4). While these anomalous efficiencies 
could be due to overestimated distances or underestimated 
moments of inertia, viewing geometries resulting in /q < 1 
could also alleviate the problem to some extent. 

We find that the fits give /q < 1 in most cases, indicating that 
we are sampling emission that is above the average over 47T . Eor 
MSPs best fit with PSPC models, /q > 1 in most cases since 
we are missing the brightest part of the emission concentrated 
at low altitudes near the polar caps. 

Eigure 7 presents the distribution of rjy for each MSP, 
excluding PSR J0610— 2100 with a suspiciously high riy ~ 13, 
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Figure 8. Ly vs. E for the 2PC MSPs with fo, estimated from the light curve 
fit with the smallest — ln(£) value. The 2PC values assuming /q = 1 are 
shown as asterisks. The E values have been derived assuming Shklovskii effect 
and Galactic acceleration corrections, as in 2PC, sometimes leading to visible 
horizontal error bars. The marker styles are as indicated in Figure 3. The dashed 
line shows 100% efficiency with the dash-dotted lines showing 10%, 1%, and 
0.1% efficiency, as labeled. The dotted line shows oc ^/e, with arbitrary 
normalization. 

(A color version of this figure is available in the online journal.) 

corrected with the /q values for the fit with the smallest — ln(£). 
These efficiencies cover a large range from ~1% to ~1200% 
with the majority clustering near 10%. 

There are three MSPs in 2PC that have greater than 100% effi- 
ciency when assuming /q=1, PSRs J0610— 2100, J0614— 3329, 
and J2302-f4442. Assuming that the distances from 2PC are not 
severely overestimated, for these MSPs, we would need /q < 
0.08, 0.47, and 0.62, respectively, to move the efficiency less 
than 100%. All of the fits for these MSPs have /q ~ 1 with esti- 
mated uncertainties that do not include the necessary /q values. 
Thus, we must conclude that if our models are correct, the high 
efficiencies of these MSPs are not strongly skewed by beaming 
effects. 

A number of MSPs now have measured masses significantly 
greater than the nominal 1.4Mq (e.g., Demorest et al. 2010; 
Ozel et al. 2012; Romani et al. 2012; Antoniadis et al. 2013). 
These larger masses could lead to larger moments of inertia 
(depending on the equation of state), increased E values, and 
lower efficiencies. 

Figure 8 plots Ly versus E as is done in 2PC but using our 
/q estimates from the light curve fitting rather than assuming 
/q=1 as was done in 2PC (asterisk points). The luminosities 
follow a roughly linear trend with E. There also seems to be 
a weak separation of TPC-preferred versus OG-preferred light 
curve fits with respect to E, with TPC-preferred for lower values 
ofE. 

8.2. Post-recycling Evolution 

Figure 9 plots the best-fit a values against the measured P. 
The Class I MSPs (filled squares and upward pointing triangles) 
display the most significant trend in a versus P, with a Pearson 
correlation coefficient of 0.5 1 and chance probability (two-sided 
p-value) of 5.6 x 10“^. The Pearson correlation coefficient tests 

A compilation of neutron star mass measurements is maintained by J. M. 
Lattimer and A. W. Steiner at http://www.stellarcollapse.org/nsmasses 
(Lattimer 2012). 
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Figure 9. Best- fit a ns. P. Only the value for the fit with the smallest — ln(£) is 
shown for each MSP. The marker styles are as indicated in Figure 3. 

(A color version of this figure is available in the online journal.) 

the linear relationship between two sets of quantities, a positive 
(negative) coefficient indicates a positive (negative) slope and 
absolute values of the coefficient closer to zero indicate weaker 
correlation. The chance probability for this trend increases to 
1.6 X 10“^ when including the Class III MSPs and 6.4 x 10“^ 
when considering all MSPs. It should be noted that the Pearson 
correlation coefficient does not incorporate the uncertainties on 
the values. In some cases, the uncertainties in a are large, so 
these trends are suggestive at best. 

A study of isolated, non-recycled pulsars by Young et al. 
(2010) suggested that as a pulsar spins down the spin and 
magnetic axes should move toward alignment (i.e., a approaches 
0°). This is in contrast to the study of Beskin & Nokhrina 
(2007), which suggested that the torque exerted by the surface 
current on the neutron star will tend to move a toward 90° 
where the energy loss is minimal. Ruderman (1991) suggested 
that spinning up a neutron star with a strong magnetic field 
(^surf ^ 10^^ G) via accretion should also move these two 
axes toward alignment. Thus, one might expect a correlation of 
increasing magnetic alignment with decreasing MSP period, 
for which we see evidence. However, the radio light curve 
and polarization data of PSR J1939-f2134, with P = 1.56 ms, 
strongly prefer an orthogonal rotator (Guillemot et al. 2012b). 
This may pose a challenge to models that suggest accretion 
should lead to nearly aligned systems or that this MSP is not a 
recycled pulsar (as has also been suggested by Ruderman 1991). 

Chen et al. (1998) noted an abundance of MSPs with both or- 
thogonal and aligned geometries compared to the non-recycled 
population, based largely on considerations of radio profiles 
more than on polarimetry or other geometrical constraints, and 
explained this as a natural consequence of the interactions be- 
tween the superfiuid neutrons and superconducting protons in- 
side a neutron star. These authors outlined two specific magnetic 
field geometries that lead to either an aligned or orthogonal ro- 
tator, and they suggested that pulsars with field configurations 
between the two cases would be somewhere between aligned 
and orthogonal. 

Our distribution of best-fit a values (Figure 3) would seem to 
be in line with the predictions of Chen et al. (1998); however. 
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Figure 10. Best- fit o; vs. E. Only the value for the ht with the smallest — ln(£) 
is shown for each MSP. The marker styles are as indicated in Figure 3. 

(A color version of this hgure is available in the online journal.) 

under their model, only the field configuration which led to 
an orthogonal rotator would be able to produce MSPs with 
the shortest periods ("^1.5 ms), and we find four MSPs with 
P < 2 ms and best-fit Ckf < 50° (though it is not clear if 
predicted polarization properties of the laSG model make it 
viable for PSR J1959-f2048 and the error on the value for 
PSR J1747— 4036 is large). Additionally, it is unclear if the 
sunspot-like geometry they assume is compatible with the VRD 
used in our simulations, even at a significant fraction of Ppc- 

Ruderman (1991) noted that a neutron star with a weak 
magnetic field spinning down should move toward an orthogonal 
alignment (in agreement with Beskin & Nokhrina 2007, though 
their findings did not depend on the magnetic field strength), 
which might explain the possible trend of increasing a with 
increasing P hinted at by the Pearson coefficient. 

Figure 10 plots the best-fit a values against the derived 
E. The Class I MSPs display the most significant trend in a 
versus logio(^), with a Pearson correlation coefficient of —0.47 
(indicating anti-correlation) and chance probability of 1.3 x 
10“^. No correlation is found when considering MSPs best 
fit by the OG and alOG models alone, but when considering 
only MSPs best fit by the TPC and alTPC models, the chance 
probability only increases to 1.5 x 10“^. When considering 
Class I and III MSPs the chance probability increases to 1.9 x 
10“^ and when considering all MSPs the chance probability is 
1.3 X 10“^ Similar to the trends in a versus P, the uncertainties 
in a are large in some cases and these trends are suggestive at 
best. 

From the Pearson coefficients, there is no clear trend in a 
as a function of E for the OG or PSPC fits in Figure 10. The 
OG model seems to prefer larger values for a, and this is to 
be expected, since there is no visible emission from this model 
when a is too low (reflecting the fact that there is no emission 
generated below the NCS). For the TPC model, however, the 
Pearson coefficient hints at a linear trend with negative slope 
between a and E. Increasing E generally points to younger age, 
and if the trend is confirmed, to a smaller a in this case. Given 
that values of ^ near 90° are usually preferred, this points to 
larger impact angles for the younger, faster- spinning MSPs in 
our sample. This corresponds to the fact that the TPC emission 
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is indeed visible for larger values of P compared to the OG 
model, since it includes low-altitude emission not present in 
the OG geometry. Furthermore, since all MSPs in this plot are 
radio-loud, we expect that the radio beams of MSPs with high 
E will generally be at higher altitudes, and therefore wider, so 
that one may probe smaller a. 

On the other hand, there may be hints of an underlying 
evolution of a with E, e.g., movement toward the equator with 
age. Is the implied, inverse trend due to unfavorable geometry, 
selection effects, or a systematic evolution of a toward larger 
values? A population synthesis approach would be necessary 
to disentangle the effects of visibility and obliquity evolution. 
We can mention in this context, though, that the relatively large 
radio beams of MSPs compared to those of the younger pulsars 
give us the opportunity to probe any supposed evolution of 
a better, though in a different magnetic field strength regime, 
given the larger parameter space over which we see both radio 
and gamma-ray emission. This is because the values of a and 
^ are much less constrained when only fitting the gamma-ray 
profiles of radio-quiet pulsars. 

8.3. MSPs with Significant Ojf-peak Emission 

An important discriminator between different emission mod- 
els is the level of emission predicted outside of the main gamma- 
ray peaks. In particular, models that restrict emission to be above 
the NCS need large \fi\to produce significant off-peak emission. 
One particular example is PSR J0613— 0200 for which Johnson 
(2011) found an OG model fit with \p\ ~ 30°. However, that 
analysis only used a hollow-cone radio model whereas we have 
included a core component for this MSP, supported by polariza- 
tion measurements, and requiring the core to be seen excludes 
this geometry. As a consequence, our best-fit OG geometry 
does not predict any off-peak emission for PSR J061 3—0200 
and misses the possible low-level peak immediately preceding 
the main gamma-ray peak. 

In 2PC, the off-peak interval was defined and analyzed for 
each pulsar with the goal of assessing the probable nature of 
any significant emission (e.g., from a pulsar wind or from 
the magnetosphere). No 2PC MSPs were found to have off- 
peak emission consistent with an origin in the pulsar wind, 
but three were found to have off-peak emission that was 
consistent with a magnetospheric origin (i.e., point-like with 
a significantly cutoff spectrum). Of these three MSPs, two 
are Class III (PSRs 10340-^4130 and J2124-3358) and one 
is Class I (PSR J2302-I-4442). The TPC model is preferred for 
PSR 12302-^442, though not significantly. The PSPC fits for 
PSRs J0340-h4130 and J2124— 3358 do predict significant off- 
peak emission, consistent with the 2PC analysis. 

There are 1 1 Class I and 4 Class II MSPs for which the best-fit 
geometries are the same to within a few degrees for the standard 
and altitude-limited TPC and OG models, neglecting laSG fits. 
Of these Class I MSPs, the likelihood prefers the TPC model for 
eight and the OG for three, but the preference is only significant 
(^15) for three MSPs (all preferring the TPC model). Of these 
Class II MSPs, the likelihood prefers the alTPC model for two 
and the alOG model for two, but never significantly. While 
some of these MSPs are faint, resulting in small log-likelihood 
difference and no clear reason for the preference, in many cases 
the preferred model is obvious from the suggestion, or not, of 
emission above the background level in the off-peak region. 
Only two of these MSPs, PSRs J0034-0534 and 12043-^1711, 
are found to have significant off-peak emission in 2PC. For both 
MSPs, the origin of the emission is unclear, while neither show 


18 


The Astrophysical Journal Supplement Series, 213:6 (54pp), 2014 July 


Johnson et al. 


any evidence for extension the significance of a cutoff in the 
off-peak spectrum is just below the threshold for classification 
as possibly magneto spheric. 

Although we have fit Class I and II MSPs with narrow 
gap models that are geometric versions of SG (TPC) and OG 
acceleration and radiation models, we find that neither provides 
the best-fit solution for all MSPs. In fact, the best-fit solutions 
divide roughly equally between the two models. The OG models 
best fit the MSP light curves that have no off-peak and high 
bridge emission, such as PSR 10030-^0451. The TPC models 
best fit the light curves having some suggestion of off-peak 
emission, such as PSR J0034— 0534. Higher off-peak emission 
may be correlated with larger ^ values that are allowed for both 
young radio-quiet pulsars and MSPs. In SG models, off-peak 
emission comes from regions below the NCS on leading-edge 
field lines. We can speculate that an emission geometry that 
would fit all MSP light curves might be an SG model in which 
the emission below the NCS either decreases toward the stellar 
surface and/or is dependent on magnetic azimuth. Such a hybrid 
model might also explain the Class I and II MSPs where the 
standard and altitude-limited TPC and OG models find similar 
geometries and the preferred fit is chosen largely based on the 
apparent presence, or not, of off-peak emission. 

8.4. Nature of the Class II MSPs 

It is still unclear what parameters give rise to the Class II 
MSPs. In particular, while the gamma-ray light curves look 
qualitatively similar to those of Class I MSPs why do the 
gamma-ray and radio components appear at nearly the same 
phase? This question has been explored by several authors 
(Johnson 2011; Venter et al. 2012; Espinoza et al. 2013), but 
no definitive conclusions have been reached. 

Harding (2005) predicted that aligned radio and HE profiles 
could occur for pulsars with P < 50 ms. With the exception 
of PSR J1823-3021A, all of the Class II MSPs have P < 
2 ms; however, PSR J 1747— 4036 has P = 1.65 ms but is 
Class I. Aside from PSRs 11939-^2134 and 11810-^1744, the 
radio profiles of Class II MSPs show essentially no linear 
polarization, in line with the prediction that caustics have low 
linear polarization (Dyks et al. 2004). Kramer et al. (1999) 
analyzed the degree of linear polarization of PSR J 1 939-^2 1 34 as 
a function of frequency and found evidence for depolarization, 
best fit by a power law with index —0.76 zb 0.07, steeper than any 
other MSP in their study. Eor PSR J1 8 10-1-1744, the situation is 
unclear as no polarimetric observations have yet been reported. 

Eigure 11 presents the E (top) and 5 lc (bottom) values 
for the 2PC MSPs separated by model class. Similar to the 
findings of Espinoza et al. (2013) but with more MSPs, a 
Kolmogorov-Smirnov test comparing the Class II MSPs (their 
A-type) to the rest of the LAT-detected MSPs in our sample 
returns a probability of 0.004 that the ^lc values are drawn 
from the same distribution and a probability of 0.02 for the E 
values. These probabilities are slightly larger than those reported 
by Espinoza et al. (2013) but still suggest that the Class II MSPs 
stand out most strongly when comparing ^lc- 

Espinoza et al. (2013) compared the radio spectral indices 
of their A-type MSPs to those of the N-type and W-type 
MSPs in their sample. A Kolmogorov-Smirnov test indicated 
a 0.6% probability that the two sets belonged to the same 
distribution. They also compared the ^lc values of their A-type 
MSPs with the rest of the gamma-ray MSP population using a 
Kolmogorov-Smirnov test and found a 0.02% probability that 
they belong to the same distribution. Espinoza et al. (2013) 



Figure 11 . Distributions of E (top) and 5 lc (bottom) by model Class. The black 
histograms filled with horizontal lines are for the Class I MSPs. The light gray 
(red in the online version) histograms filled with downward diagonal lines from 
left to right are for the Class II MSPs. The dark gray filled histograms (blue in 
the online version, filled with upward diagonal lines from left to right) are for 
the Class III MSPs. 

(A color version of this figure is available in the online journal.) 

also pointed out that three of the A-type MSPs are known to 
emit giant radio pulses, as is PSR J0218-f4232, a W-type MSP. 
Significant gamma-ray pulsations have recently been detected 
from PSR B 1821— 24 in the globular cluster M28 (Wu et al. 
2013; Johnson et al. 2013), which is also known to emit giant 
pulses. The first gamma-ray, and X-ray, light curve peaks occur 
at the same phase from which the giant pulses have been detected 
but the identification of this MSP as A-type is unclear. Giant 
pulses have also been detected from the Crab pulsar (Heiles 
et al. 1970), the only non-recycled pulsar with the main radio 
and gamma-ray peaks occurring at the same phase. 

It is not clear what physical characteristics divide MSPs into 
the three model classes we have identified. Class I and II do not 
show a preference for any particular viewing geometry, although 
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Pulse Phase 

Figure 12. Best- fit gamma-ray (top) and radio (bottom) light curves for PSR J0023+0923. The dashed line (solid green in the online version) is the OG model, and 
best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and best-ht parameters are given in Table 6. In both panels, 
the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 


Class III MSPs all have ^ > a. Since very few Class I and II 
MSPs were expected originally, perhaps class membership is 
governed by other characteristics that determine whether narrow 
gaps can form in these MSPs, such as surface magnetic field 
structure. 

The laSG model is usually not preferred. It requires small gap 
widths, with the gamma-ray widths generally smaller than the 
radio widths. The preferred positions of the peak emission seem 
to be about 2 7?ns with the best-fit a and ^ generally exceeding 
40° and 70°, respectively. 

The large number of gamma-ray MSPs discovered by the 
Fermi LAT has raised some intriguing questions not only about 
MSP emission and evolution, but also about pulsar emission 
physics in general. Continued study of these sources is sure to 
yield more surprises and lead to more answers. 
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APPENDIX A 
ALL LIGHT CURVE EITS 

Eigure 12 presents the observed and best-fit light curves for 
PSR 10023-^0923. This is a 3.05 ms pulsar in a binary system 
with a 0.14 day orbital period discovered in radio observations 
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Figure 13. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J0030+0451. The dashed line (solid green in the online version) is the OG model, and 
best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and best-ht parameters are given in Table 6. In both panels, 
the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 14. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR J0034— 0534. The dashed line (solid green in the online version) is the alOG model, 
and best-ht parameters are given in Table 9. The solid gray line (pink in the online version) is the alTPC model, and best-ht parameters are given in Table 8. The 
dash-dotted, light gray line (solid red in the online version) is the laSG model, and best-ht parameters are given in Table 10. In both panels, the dotted horizontal lines 
are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 


21 


The Astrophysical Journal Supplement Series, 213:6 (54pp), 2014 July 


Johnson et al. 



Figure 15. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR JOlOl— 6422. The dashed line (solid green in the online version) is the OG model, and 
best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and best-ht parameters are given in Table 6. In both panels, 
the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 16. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR JOl 02+4839. The dashed line (solid green in the online version) is the OG model, and 
best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and best-ht parameters are given in Table 6. In both panels, 
the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 17. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J02 18+4232. The dashed line (solid green in the online version) is the OG model, and 
best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and best-ht parameters are given in Table 6. In both panels, 
the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 18. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR J0340+4130. The solid gray line (blue in the online version) is the PSPC model, and 
best-ht parameters are given in Table 1 1. In both panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 19. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J0437— 4715. The dashed line (solid green in the online version) is the OG model, and 
the best-fit parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-fit parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this figure is available in the online journal.) 




Figure 20. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J0610— 2100. The dashed line (solid green in the online version) is the OG model, and 
the best-fit parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-fit parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this figure is available in the online journal.) 
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Figure 21. Best- fit gamma-ray (top) and radio (bottom) light curves for PSR J0613— 0200. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best- fit parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this figure is available in the online journal.) 



Figure 22. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J0614— 3329. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 23. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J0751+1807. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 24. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR J1 024—07 19. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. The 
best-ht geometries are the same resulting in the same radio, model light curve resulting in the OG ht being drawn on top of the TPC ht. In both panels, the dotted 
horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 25. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1 124—3653. The dashed line (solid green in the online version) is the OG model, and 
the best- fit parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. The 
best-ht geometries are the same resulting in the same radio, model light curve resulting in the OG fit being drawn on top of the TPC fit. In both panels, the dotted 
horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 


of unassociated LAT sources with pulsar-like characteristics 
(J. W. T. Hessels et ah, in preparation). This MSP was first 
announced as a gamma-ray pulsar in 2PC. Both the OG and 
TPC models fit the gamma-ray light curve well, though neither 
best-fit geometry matches the smaller radio peak near phase 0.4. 
Breton et al. (2013) modeled the optical light curve of this MSP 
and found i = 58° zb 14° and were able to exclude i < 40°. 
This does agree with the value of ^ for our TPC fit, within 
uncertainties, and it better matches the lack of observed radio 
eclipses. Only the TPC confidence contours can accommodate 
^ = 58°, with a slightly lower value of a than that in Table 6. 

Figure 13 presents the observed and best-fit light curves for 
PSR J0030-f0451. This is an isolated pulsar with a 4.87 ms 
spin period first discovered by Lommen et al. (2000). Gamma- 
ray pulsations from this MSP were first reported by Abdo 
et al. (2009c). Lommen et al. (2000) also attempted to fit the 
polarimetric data for this MSP but could not uniquely constrain 
the viewing geometry. They did suggest that the pulsar was 
either nearly aligned or nearly orthogonal, and Bogdanov et al. 
(2008) ruled out the nearly aligned geometry by modeling the 
thermal X-ray pulsations. Fits with the TPC and OG model 
both find values of a near 90°, but the OG fit is significantly 
preferred by the likelihood as the TPC model overpredicts the 
level of off-peak emission and only the OG model produces a 
radio interpulse. 

Figure 14 presents the observed and best-fit light curves of 
PSR J0034— 0534. This is a 1.88 ms pulsar in a binary system 
with a 1.6 day orbital period discovered by Bailes et al. (1994). 
Gamma-ray pulsations from this MSP were first reported by 
Abdo et al. (2010b). The radio emission from this MSP shows 


very low levels of linear polarization, consistent with being 
unpolarized (Stairs et al. 1999). The alTPC, alOG, and laSG 
models all fit the observed light curves well, but only the 
laSG and alTPC models predict significant levels of off-peak 
emission. The laSG model does the best at matching the relative 
radio peak heights but underestimates the gamma-ray peak at 
phase 0.8. 

Figure 15 presents the observed and best-fit light curves 
of PSR JOlOl— 6422. This is a 2.57 ms pulsar in a binary 
system with a 1.8 day orbital period discovered by Kerr et al. 
(2012), who also announced this as a gamma-ray MSP, in 
radio observations of unassociated LAT sources with pulsar- 
like characteristics. Neither the TPC or OG model is able to 
simultaneously match the gamma-ray and radio light curves. 
The OG model fits the gamma-rays reasonably well and is the 
only model to predict more than one radio peak. 

Figure 16 presents the observed and best-fit light curves of 
PSR J0102-f4839. This is a 2.96 ms pulsar in a binary system 
with a 1.7 day orbital period discovered in radio observations 
of unassociated LAT sources with pulsar-like characteristics 
(P. Bangale et al., in preparation). This MSP was first announced 
as a gamma-ray pulsar in 2PC. Both the OG and TPC model 
are able to reproduce the observed gamma-ray light curve well, 
but both also predict a radio lag that is too small. The best-fit 
geometry corresponding to the TPC fit does match the observed 
radio profile better. 

Figure 17 presents the observed and best-fit light curves of 
PSR J0218-F4232. This is a 2.32 ms pulsar in a binary system 
with a 2.0 day orbital period discovered by Navarro et al. (1995). 
Low- significance pulsations (^3a) using EGRET data were 
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Figure 26. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1 125— 5825. The dashed line (solid green in the online version) is the OG model, and 
the best-fit parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-fit parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this figure is available in the online journal.) 



Figure 27. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1231 — 141 1. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best- fit parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 28. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1446— 4701. The dashed line (solid green in the online version) is the OG model, and 
the best-fit parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-fit parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this figure is available in the online journal.) 


reported by Kuiper et al. (2000), with a firm pulsed detection 
using LAT data first reported by Abdo et al. (2009a). Stairs 
et al. (1999) reported polarization measurements for this MSP 
at 410 and 610 MHz and found o' = 8° di 11°(15°) for the 410 
(610) MHz data, with unconstrained from RVM fits. Neither 
the TPC or OG model fits exactly match the gamma-ray light 
curve of this MSP, but only the TPC model fit matches the radio 
reasonably well. The best-fit a from the TPC model agrees 
with the findings of Stairs et al. (1999) within statistical and 
systematic uncertainties, while that from the OG model does 
not, largely due to the fact that only accepting emission from 
above the NCS tends to exclude low-a solutions. 

Figure 18 presents the observed and best-fit light curves of 
PSR J0340-f4130. This is an isolated 3.30 ms pulsar discovered 
in radio observations of unassociated LAT sources with pulsar- 
like characteristics (P. Bangale et al., in preparation). This MSP 
was first announced as a gamma-ray pulsar in 2PC. The PSPC 
model is able to reproduce both the radio and gamma-ray light 
curves of this MSP reasonably well. The gamma-ray fit does 
not find two well- separated peaks, nor is the main peak sharp 
enough. 

Figure 19 presents the observed and best-fit light curves of 
PSR J0437— 4715. This is a 5.76 ms pulsar in a binary system 
with a 5.7 day orbital period first discovered by Johnston et al. 
(1993). Gamma-ray pulsations were first reported by Abdo 
et al. (2009a). Manchester & Johnston (1995) and Navarro 
et al. (1997) both presented polarimetric observations of this 
MSP, including circular polarization with sense reversal near 
the profile center indicative of a core component. Thus, in our 


fits we used both a hollow-cone and core component and require 
that the best-fit ^ intersects the core component. Manchester & 
Johnston (1995) attempted to constrain the viewing geometry 
of PSR J0437— 4715 using RVM fits but were unsuccessful. 
They did present RVM curves for a = 145° and ^ = 140° 
which seemed to represent the data well (these values would 
correspond to o' = 35° and ^ = 40° in our fits). Gil & 
Krawczyk (1997) applied the relativistic RVM of Blaskiewicz 
et al. (1991) to the same data and found a = 20° and ^ = —4°. 
These authors first found geometries which reproduced the 
basic profile characteristics and then calculated the predicted 
polarization properties and compared to the observations, so 
the values they report are not really fits and do not come with 
any uncertainty estimates. Bogdanov et al. (2007) modeled the 
thermal X-ray light curves of this MSP and found an offset- 
dipole geometry was necessary to match the data. They assumed 
that the spin and orbit axis were aligned, and thus used ^ 
= 42°, from the timing analyses of van Straten et al. (2001) 
and Hotan et al. (2006) who found orbital inclination angles 
i = 42?75 zb 0?09 and 42?5 zb 0?2, respectively. With the 
offset dipole and this assumed ^ Bogdanov et al. (2007) found 
acceptable fits to the observed X-ray light curves. This value of 
^ would argue against the geometry found by Gil & Krawczyk 
(1997). The predicted radio light curves from both the TPC and 
OG fits match the observed profile well (with this binning many 
of the lower-level features of the profile are smoothed out), but 
neither exactly matches the gamma-ray light curve. The OG 
model best fit finds ^ = 46°, in better agreement with the orbital 
inclination measurements than the value of 64° found by the 
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Figure 29. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1 5 14— 4946. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best- fit parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 



Pulse Phase 


Figure 30. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1 600— 3053. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 31. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1614— 2230. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best- fit parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 



Pulse Phase 


Figure 32. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1658— 5324. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 33. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J17 13+0747. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 



Figure 34. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR J1741+1351. The solid gray line (blue in the online version) is the PSPC model, and the 
best-ht parameters are given in Table 1 1. In both panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 35. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1744— 1134. The solid gray line (blue in the online version) is the PSPC model, and the 
best-ht parameters are given in Table 1 1. In both panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 36. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR J1747— 4036. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. The 
best-ht geometries are the same resulting in the same radio, model light curve resulting in the OG ht being drawn on top of the TPC ht. In both panels, the dotted 
horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 37. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1810+1744. The dashed line (solid green in the online version) is the alOG model, and 
the best-fit parameters are given in Table 9. The solid gray line (pink in the online version) is the alTPC model, and the best-fit parameters are given in Table 8. The 
dash-dotted, light gray line (solid red in the online version) is the laSG model, and the best- fit parameters are given in Table 10. In both panels, the dotted horizontal 
lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 38. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR J1 823— 3021 A. The dashed line (solid green in the online version) is the alOG model, and 
the best-fit parameters are given in Table 9. The solid gray line (pink in the online version) is the alTPC model, and the best-fit parameters are given in Table 8. The 
dash-dotted, light gray line (solid red in the online version) is the laSG model, and the best-ht parameters are given in Table 10. In both panels, the dotted horizontal 
lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 39. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1 858— 2216. The solid gray line (blue in the online version) is the PSPC model, and the 
best-ht parameters are given in Table 1 1. In both panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 40. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR J1902— 5105. The dashed line (solid green in the online version) is the alOG model, and 
the best-ht parameters are given in Table 9. The solid gray line (pink in the online version) is the alTPC model, and the best-ht parameters are given in Table 8. The 
dash-dotted, light gray line (solid red in the online version) is the laSG model, and the best-ht parameters are given in Table 10. In the radio panel, the altitude-limited 
models predict the same radio light curve resulting in the alOG ht being drawn on top of the alTPC ht. In both panels, the dotted horizontal lines are the estimated 
background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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TPC fit. In the future, fits to the radio and gamma-ray light 
curve of this MSP with an offset-dipole geometry will be useful 
to more directly compare with the results of Bogdanov et al. 
(2007). 

Figure 20 presents the observed and best-fit light curves of 
PSR J0610— 2100. This is a 3.86 ms pulsar in a binary system 
with a 0.29 day orbital period discovered by Burgay et al. (2006). 
Gamma-ray pulsations from this MSP were first reported by 
Espinoza et al. (2013). Both the OG and TPC models are able 
to reproduce the gamma-ray light curve of this MSP reasonably 
well, and both predict that the light curve should consist of 
two peaks, though more statistics are necessary to confirm the 
existence of a second peak. The predicted radio light curve from 
the TPC fit matches the main radio peak better, but both models 
predict two radio peaks, with different relative heights, though 
the second peak is not at the correct phase in either model. 

Figure 21 presents the observed and best-fit light curves 
of PSR J0613— 0200. This is a 3.06 ms pulsar in a binary 
system with a 1.2 day orbital period discovered by Lorimer 
et al. (1995). Gamma-ray pulsations from this MSP were first 
reported by Abdo et al. (2009a). Polarimetric observations for 
this MSP have been reported by Xilouris et al. (1998), Stairs 
et al. (1999), Ord et al. (2004), Manchester & Han (2004), 
and Yan et al. (2011) with all but Manchester & Han (2004) 
reporting circular polarization with sense reversal through the 
main pulse component. As such, we use a hollow-cone and core 
radio beam in our fits and require that the best-fit ^ intersects 
the core component. Both the OG and TPC model reproduce the 
gamma-ray light curve well, but only the TPC model matches 
the possible precursor to the main peak at phase 0.2. Neither 
best-fit geometry leads to a wide enough radio pulse or to the 
correct peak multiplicity. 

Figure 22 presents the observed and best-fit light curves of 
PSR J0614— 3329. This is a 3.15 ms pulsar in a binary system 
with a 54 day orbital period discovered by Ransom et al. (201 1), 
who also reported the first detection of gamma-ray pulsations 
from this MSP, in radio observations of FAT unassociated 
sources with pulsar-like characteristics. While both the OG and 
TPC model predict a two-peaked gamma-ray light curve, neither 
quite matches the observed peak separation nor the peak-height 
ratio. Both best-fit geometries do predict two radio peaks; but 
the smaller peak is not at the correct phase, and the main peak 
is not wide enough, though the TPC fit does come close. 

Figure 23 presents the observed and best-fit light curves 
for PSR J0751-hl807. This is a 3.48 ms pulsar in a binary 
system with a 0.26 day orbital period discovered by Fundgren 
et al. (1995). Gamma-ray pulsations from this MSP were first 
reported by Abdo et al. (2009a). Xilouris et al. (1998) reported 
polarimetric observations of this MSP but make no predictions 
on the viewing geometry. Neither the TPC nor OG model 
matches the three peaks evident in the gamma-ray light curve 
and only the best-fit geometry corresponding to the OG fit 
predicts two closely spaced radio peaks. 

Figure 24 presents the observed and best-fit light curves of 
PSR J1024— 0719. This is an isolated 5.16 ms pulsar discovered 
by Bailes et al. (1997). Gamma-ray pulsations from this MSP 
were first reported by Espinoza et al. (2013). Yan et al. (2011) 
reported polarimetric observations of this MSP, the position- 
angle swing of the linear polarization is fiat, suggesting our line 
of sight clips the edge of a cone. They also noted weak circular 
polarization with sense-reversal through the main component. 
The circular polarization appears to have substructure and at 
least three overlapping components, so identification of a core 
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beam is uncertain. Since Yan et al. (2011) do not specifically 
call this out as evidence for a core component, the results we 
report in Section 7 only use a hollow-cone beam for the radio 
model. However, when we include a core beam and require that 
the best-fit ^ intersects this component, the likelihood finds the 
same best-fit geometry. The best-fit model does predict that the 
core component is much weaker than the cone. Neither the OG 
nor the TPC model reproduces the observed gamma-ray light 
curve well with both predicting two peaks. However, this MSP 
is relatively faint in gamma rays and the estimated background 
level may be too high. Future modeling with more statistics will 
be interesting. 

Figure 25 presents the observed and best-fit light curves of 
PSR J1 124— 3653. This is a 2.41 ms pulsar in a binary system 
with a 0.23 day orbital period discovered in radio observations 
of unassociated FAT sources with pulsar-like characteristics 
(P. Bangale et al., in preparation). This MSP was first announced 
as a gamma-ray pulsar in 2PC. Both the TPC and OG models 
reproduce the gamma-ray light curve well but only match the 
main radio peak. 

Figure 26 presents the observed and best-fit light curves of 
PSR J1 125— 5825. This is a 3.10 ms pulsar in a binary system 
with a 76 day orbital period discovered by Bates et al. (2011). 
Keith et al. (2012) reported a detection of pulsed gamma rays 
from this MSP with a significance of 4.9a. They also reported 
polarimetric observations of this MSP and used an RVM fit 
of only the main pulse component to predict a = 128° and 
P = —21° with no quoted uncertainties. Both the TPC and OG 
model can qualitatively reproduce the observed gamma-ray light 
curve, though more statistics would lead to a more well-defined 
profile. Both models match the main radio peak but only the 
OG model fit geometry predicts more than 1 peak, though the 
second peak is not at the correct phase. The viewing geometry 
found by Keith et al. (2012) would correspond to o' = 52° and 
^ = 73° in our fits. Both the TPC and OG models find values of 
^ which are compatible with their estimate, but only the TPC a 
value is consistent with 52° within the quoted uncertainties. 

Figure 27 presents the observed and best-fit light curves of 
PSR J1 23 1 — 1411. This is a 3.68 ms pulsar in a binary system 
with a 1 .9 day orbital period discovered by Ransom et al. (201 1), 
who also reported the first detection of gamma-ray pulsations 
from this MSP, in radio observations of unassociated FAT 
sources with pulsar-like characteristics. Both the TPC and OG 
models can reproduce approximately correct phase separation 
for the gamma-ray peaks, with the OG model matching the 
shapes better. The best-fit geometry corresponding to the TPC 
fit matches the radio profile better than the OG model fit, which 
predicts a nearly orthogonal geometry leading to two radio 
peaks, one from each pole. 

Figure 28 presents the observed and best-fit light curves for 
PSR J1446— 4701. This is a 2.19 ms pulsar in a binary system 
with a 0.28 day orbital period discovered by Keith et al. (2012), 
who also reported the first detection of gamma-ray pulsations 
from this MSP. Both the TPC and OG model reproduce the 
gamma-ray light curve of PSR J1446— 4701 well, with the TPC 
model matching the position of the peak slightly better. Both 
models also match the radio profile well, though with slightly 
larger predicted phase lag than what is observed. 

Figure 29 presents the observed and best-fit light curves of 
PSR J1 5 14—4946. This is a 3.58 ms pulsar in a binary system 
with a 1 .9 day orbital period discovered in a radio observations 
of unassociated FAT sources with pulsar-like characteristics 
and first announced by Kerr et al. (2012) (with more details in 
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Figure 41. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1939+2134. The dashed line (solid green in the online version) is the alOG model, and 
the best-fit parameters are given in Table 9. The solid gray line (pink in the online version) is the alTPC model, and the best-fit parameters are given in Table 8. The 
dash-dotted, light gray line (solid red in the online version) is the laSG model, and the best- fit parameters are given in Table 10. In both panels, the dotted horizontal 
lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 42. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J1959+2048. The dashed line (solid green in the online version) is the alOG model, and 
the best-fit parameters are given in Table 9. The solid gray line (pink in the online version) is the alTPC model, and the best-fit parameters are given in Table 8. The 
dash-dotted, light gray line (solid red in the online version) is the laSG model, and the best-ht parameters are given in Table 10. In both panels, the dotted horizontal 
lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 


37 


The Astrophysical Journal Supplement Series, 213:6 (54pp), 2014 July 


Johnson et al. 



Figure 43. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J20 17+0603. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 



Figure 44. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR J2043+1711. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 45. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J2047+1053. The dashed line (solid green in the online version) is the OG model, and the 
best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best- fit parameters are given in Table 6. The best-fit 
geometries are the same resulting in the same radio, model light curve resulting in the OG fit being drawn on top of the TPC fit. In both panels, the dotted horizontal 
lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 



Figure 46. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J205 1—0827. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 47. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J2 124— 3358. The solid gray line (blue in the online version) is the PSPC model, and the 
best-ht parameters are given in Table 1 1. In both panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 48. Best-ht gamma-ray (top) and radio (bottom) light curves for PSR J2214+3000. The solid gray line (blue in the online version) is the PSPC model, and the 
best-ht parameters are given in Table 1 1. In both panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 49. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J22 15+5 135. The dashed line (solid green in the online version) is the OG model, and the 
best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-fit parameters are given in Table 6. The best-fit 
geometries are the same resulting in the same radio, model light curve resulting in the OG fit being drawn on top of the TPC fit. In both panels, the dotted horizontal 
lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Figure 50. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J2241— 5236. The dashed line (solid green in the online version) is the OG model, and 
the best- fit parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. The 
best-ht geometries are nearly the same resulting in similar radio, model light curve resulting in the OG ht being drawn on top of the TPC ht. In both panels, the dotted 
horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 
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Pulse Phase 

Figure 51. Best-fit gamma-ray (top) and radio (bottom) light curves for PSR J2302+4442. The dashed line (solid green in the online version) is the OG model, and 
the best-ht parameters are given in Table 7. The solid gray line (pink in the online version) is the TPC model, and the best-ht parameters are given in Table 6. In both 
panels, the dotted horizontal lines are the estimated background levels given in Table 5. 

(A color version of this hgure is available in the online journal.) 


F. Camilo et al., in preparation). This MSP was first announced 
as a gamma-ray pulsar in 2PC. Both the TPC and OG models 
match the gamma-ray light curves of this MSP well, but neither 
model is able to reproduce the two closely spaced radio peaks 
but the qualitative properties of the predicted radio emission 
agree well with observations. 

Figure 30 presents the observed and best-fit light curves of 
PSR J1600— 3053. This is a 3.60 ms pulsar in a binary system 
with a 14 day orbital period discovered by Jacoby et al. (2007). 
Ord et al. (2004) and Yan et al. (2011) reported polarimetric 
observations for this MSP and noted sense reversal in the circular 
polarization. As such, we have used both a hollow-cone beam 
and a core beam when modeling the light curves of this pulsar 
and required the best-fit ^ to intersects the core component. 
The best-fit geometries corresponding to fits with the TPC and 
OG models both reproduce the observed radio profile well. 
However, neither fit is able to properly match the first gamma- 
ray peak. This MSP is relatively faint in gamma rays and more 
statistics may lead to a sharper profile and drive the likelihood 
to a geometry which better matches the data. 

Figure 3 1 presents the observed and best-fit light curves of 
PSR J1614— 2230. This is a 3.15 ms pulsar in a binary system 
with a 8.7 day orbital period discovered by Crawford et al. 
(2006) in radio observations of unassociated EGRET sources. 
Gamma-ray pulsations from this MSP were first reported by 
Abdo et al. (2009a). Demorest et al. (2010) used the Shapiro 
delay to measure the mass of this pulsar to be (1.97 zb 0.04) Mq 
and the inclination angle with respect to the orbit axis to be 
i = 88 ? 17 zb 0?02. If we assume that the spin and orbit axes 
are aligned, this can be taken as a constraint on ^ , but we only 


compare our best-fit values to this and do not use it to limit 
the likelihood analysis. Both the OG and TPC model match 
the observed gamma-ray light curve well, but neither exactly 
matches the radio profile. While both best-fit geometries predict 
radio emission from two poles, the relative peak heights and 
phases are wrong. Both models find values of ^ near 90° with 
the OG model agreeing best with ^ = ( 88 ^. 5 )°. 

Figure 32 presents the observed and best-fit light curves of 
PSR J1658— 5324. This is an isolated 2.44 ms pulsar discovered 
in radio observations of unassociated LAT sources with pulsar 
like characteristics, first announced by Kerr et al. (2012; with 
more details in F. Camilo et al., in preparation). This MSP 
was first announced as a gamma-ray pulsar in 2PC. Both 
the OG and TPC model match the main gamma-ray peak of 
PSR J1658— 5324 well, with the TPC model predicting two 
very closely spaced peaks. The best-fit geometries for both fits 
are able to match the highest radio peak, but they miss the broad 
peak preceding this and the small peak near phase 0 . 8 . 

Figure 33 presents the observed and best-fit light curves 
of PSR J1713-f0747. This is a 4.57 ms pulsar in a binary 
system with a 68 day orbital period discovered by Foster 
et al. (1993). Gamma-ray pulsations from this MSP were first 
reported by Espinoza et al. (2013). Several authors have reported 
polarimetric observations of PSR J1713-f0747 (e.g., Xilouris 
et al. 1998; Stairs et al. 1999; Ord et al. 2004; Yan et al. 
2011 ) and all measure some degree of circular polarization, 
with sense reversal in some cases, but none have attempted 
RVM fits for this MSP. We use both a hollow-cone and core 
beam to model the radio light curve of this MSP and require 
that the best-fit ^ intersects the core beam. This pulsar is 
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one of the best-timed and several authors have used Shapiro 
delay measurements to constrain the orbital inclination angle. 
Camilo et al. (1994) were the first to report such measurements 
for this MSP, and they estimated i = 70°, though the fits 
were not well constrained. Splaver et al. (2005) and Hotan 
et al. (2006) reported more precise timing measurements and 
found i = 71?9 zb 1?8 and (76. 6 ^. 20 )°’ respectively. As with 
PSR J1 6 14— 2230, these values can be used as constraints on 

assuming alignment of the spin and orbit axes, but we only 
compare our best-fit values to these and do not include them as 
limits in the likelihood fit. Both the TPC and OG model are able 
to reproduce the observed gamma-ray light curve of this MSP 
well. The TPC model matches the position of the main peak 
well with the OG model predicting two closely spaced peaks, 
but more statistics will be necessary to determine if this MSP 
has one or two gamma-ray peaks. The best-fit geometries for 
both fits reproduce the observed radio profile well. Our TPC 
fit finds ^ = 68° zb 4° and our OG fit finds (65^.4^)°. Both 
fits agree well, within uncertainties, with the constraints from 
timing. 

Figure 34 presents the observed and best-fit light curves of 
PSR J1741-hl351. This is a 3.75 ms pulsar in a binary system 
with a 16 day orbital period discovered by Jacoby et al. (2007). 
Gamma-ray pulsations from this MSP were first reported by 
Espinoza et al. (2013). The PSPC model is able to reproduce the 
single gamma-ray peak of this MSP well, though the predicted 
peak is not sharp enough and it is unclear if the predicted 
off-peak emission is seen with the LAT. The best-fit geometry 
matches the main radio peak well but does not reproduce the 
smaller peak at phase 0.9. 

Figure 35 presents the observed and best-fit light curves 
of PSR J1744— 1134. This is an isolated 4.07 ms pulsar 
discovered by Bailes et al. (1997). Gamma-ray pulsations from 
this MSP were first reported by Abdo et al. (2009a). Polarimetric 
observations of this pulsar all present fiat position- angle swings 
favoring a cone beam (Xilouris et al. 1998; Stairs et al. 1999; 
Ord et al. 2004; Yan et al. 201 1). Additionally, Yan et al. (201 1) 
suggest that > 0 based on the assumption that the position- 
angle swing in the main component is a continuation of that in 
the precursor. The PSPC model is able to reproduce the observed 
gamma-ray light curve well, and the best-fit geometry matches 
the main radio peak but not the low-level precursor. 

Figure 36 presents the observed and best-fit light curves of 
PSR J1747— 4036. This is an isolated 1.65 ms discovered in 
radio observations of unassociated FAT sources with pulsar- 
like spectra, first announced by Kerr et al. (2012; with more 
details presented in F. Camilo et al., in preparation). This MSP 
was first announced as a gamma-ray pulsar in 2PC. The gamma- 
ray peak and main radio peak are nearly aligned in phase but 
there is no gamma-ray component matching the second radio 
peak and the polarization properties of this MSP (see F. Camilo 
et al., in preparation) are not in line with most other Class II 
MSPs. Both the OG and TPC models reproduce the gamma-ray 
light curve well and match the main radio peak, but neither can 
reproduce the second radio peak. 

Figure 37 presents the observed and best-fit light curves of 
PSR J1 8 10-^1744. This is a 1.66 ms pulsar in a binary system 
with a 0.15 day orbital period discovered in radio observations 
of unassociated FAT sources with pulsar-like characteristics 
(J. W. T. Hessels et al., in preparation). This MSP was first 
announced as a gamma-ray pulsar in 2PC. The broad gamma- 
ray and radio peaks occur at the same phase, but the shapes do 
not match exactly and it is unclear if the smaller radio peak at 
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phase 0.9 has a matching feature in gamma rays. The alTPC, 
alOG, and laSG models all match the qualitative features of 
the gamma-ray and radio light curves, with the alOG model 
coming the closest to reproducing the three gamma-ray peaks 
in the observed light curve and only the alTPC model matching 
the smaller radio peak. The optical light curve of this MSP 
has been modeled by Breton et al. (2013) and Schroeder & 
Halpern (2014). The first authors found / = 48° zb 7°, though 
they note that their fit was poor and unable to match all of the 
data points. The latter authors used two different models, finding 
i = 56?75 zb 2?25 and 54?75 zb 2.75, compatible with the results 
of Breton et al. (2013); however, they were unable to produce 
fits that led to a realistic neutron star mass or X-ray efficiency. 
This does not agree well with any of our fit values of ^ for this 
MSP and only the confidence contours for the laSG fit allow for 
solutions with ^ 55°. 

Figure 38 presents the observed and best-fit light curves 
of PSR J1823— 3021A. This is an isolated 5.44 ms pulsar 
discovered in the globular cluster NGC 6624 by Biggs et al. 
(1994). Gamma-ray pulsations from this MSP were first reported 
by Freire et al. (2011) who suggested, based on the inferred 
gamma-ray efficiency, that the P (and thus E\ see Table 4) 
of this MSP was largely intrinsic and, therefore, this is an 
extremely luminous gamma-ray pulsar. Stairs et al. (1999) 
reported polarimetric observations of this MSP, which show no 
position-angle swing across the pulse and, in fact, are consistent 
with a mean polarization of 0%. The estimated background level 
in the top panel of Figure 38 may be overestimated, making 
it unclear how large the off-peak interval is and how many 
components the gamma-ray light curve has. The alTPC, alOG, 
and laSG models all match the main gamma-ray peak well but 
predict different peak multiplicities. The predicted radio profiles 
suffer from the same issue. 

Figure 39 presents the observed and best-fit light curves of 
PSR J1 858—2216. This is 2.38 ms pulsar in a binary system 
with a 46 day orbital period discovered in radio observations of 
unassociated FAT sources with pulsar-like characteristics 
(S. Sanpa-Arsa et al., in preparation). This MSP was first an- 
nounced as a gamma-ray pulsar in 2PC. The PSPC model re- 
produces the gamma-ray peak well, though the predicted peak 
is not sharp enough and the predicted level of off-peak emission 
may not be seen with FAT. The best-fit geometry matches the 
main radio peak well but does not match the possible low-level 
feature between phases of 0.2 and 0.4, though it is unclear from 
our background level estimate whether or not this feature is real. 

Figure 40 presents the observed and best-fit light curves of 
PSR J1902— 5105. This is a 1.74 ms pulsar in a binary system 
with a 2.0 day orbital period discovered in radio observations 
of unassociated FAT sources with pulsar-like characteristics, 
first announced by Kerr et al. (2012; with more details in 
F. Camilo et al., in preparation). This MSP was first announced 
as a gamma-ray pulsar in 2PC. The radio emission from this 
pulsar shows little to no polarization (F. Camilo et al., in 
preparation). The alOG, alTPC, and laSG models all match 
the observed gamma-ray and radio light curves of this MSP 
well. None of the models match the sharpness of the gamma- 
ray peak near phase 0.9 properly, but this may be due to the 
likelihood putting some emphasis at matching the apparent 
bridge emission. None of the models exactly match the shapes 
of the radio peaks but the alTPC and alOG models do better at 
predicting the proper peak-height ratio. 

Figure 41 presents the observed and best-fit light curves of 
PSR J1939-F2134 (B1937-h21, the first MSP ever discovered. 
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Backer et al. 1982). This is an isolated 1.56 ms pulsar. 
Gamma-ray pulsations from this MSP were first reported by 
Guillemot et al. (2012b) who also predicted the viewing ge- 
ometry from RVM fits to be o' = 89° and = —3° with a 
radio emission altitude of 0.65 Rpc- Of the Class II MSPs for 
which polarimetric observations have been reported, this is the 
only one with significant levels of polarized emission. Several 
authors have reported polarimetric observations for this MSP 
(Ashworth et al. 1983; Stinebring 1983; Stinebring & Cordes 
1983; Thorsett & Stinebring 1990; Xilouris et al. 1998; Stairs 
et al. 1999; Ord et al. 2004; Yan et al. 2011), and all agree on the 
basic properties of the emission. The near 180° phase separation 
of the peaks suggests an orthogonal rotator viewed near the spin 
equator (i.e., both a and ^ near 90°), Stinebring (1983) noted that 
the polarization is consistent with this geometry if the depolar- 
ization near the inner peak edges was a sign of orthogonal mode 
switches. This depolarization for an orthogonal rotator would 
also be expected from outer-magnetospheric emission models 
(Dyks et al. 2004). Thorsett & Stinebring (1990) noted that the 
narrowness of the peaks is at odds with the typical interpretation 
of a core beam for an orthogonal rotator. A caustic origin of the 
emission does naturally explain the narrow peaks, but so would 
partially filled cone beams (Kramer et al. 1999). The alOG, 
alTPC, and laSG models all predict the observed gamma-ray 
and radio light curves well, though none exactly matches the 
phase of the gamma-ray peak near phase 0.5. Only the best-fit 
geometry corresponding to the alTPC fit agrees well with the 
RVM fits of Guillemot et al. (2012b), though the confidence 
contours for the alOG fit are not simply connected and do allow 
a fit consistent with their values within the 95% confidence-level 
region. 

Figure 42 presents the observed and best-fit light curves of 
PSR J1959-F2048 (B 1957-^20, the first “black- widow” pulsar 
ever discovered, Fruchter et al. 1988). This is a 1.61 ms pulsar 
in a binary system with a 0.38 day orbital period and a very 
low-mass (0.022 Mq) companion. The companion is thought 
to have reached such a low mass from ablation by the pulsar 
wind, thus the “black widow” moniker. Gamma-ray pulsations 
were first reported from PSR J1959-I-2048 by Guillemot et al. 
(2012b). Thorsett & Stinebring (1990) reported polarimetric 
observations of this MSP, finding very-low levels of linear 
polarization (<2% of the total intensity). Reynolds et al. (2007) 
used optical observations of the companion of PSR 11959-^2048 
to constrain the inclination of the orbit to be i = 65° ± 2°, 
assuming a pulsar mass between 1.3 and 1.9 Mq. The alOG, 
alTPC, and laSG models are all able to reproduce the observed 
gamma-ray light curves well and match the first and third radio 
peaks, but none of the models can match the second radio peak 
which does not have a corresponding feature in the gamma- 
ray light curve. The second radio peak, near phase 0.9, is not 
present at lower frequencies. None of the models find best-fit ^ 
values that agree well with the orbital inclination of Reynolds 
et al. (2007), though the laSG value is consistent within our 
estimated uncertainties. 

Figure 43 presents the observed and best-fit light curves for 
PSR J2017-F0603. This is 2.90 ms pulsar in a binary system 
with a 2.2 day orbital period discovered by Cognard et al. 
(2011), who also reported the first detection of gamma-ray 
pulsations from this MSP, in radio observations of unasso- 
ciated LAT sources with pulsar-like characteristics. The OG 
and TPC models both reproduce the gamma-ray peaks well, 
though the TPC model over predicts the off-peak emission and 
predicts a precursor at phase 0.3, which is not present in the 
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data. Neither model results in a best-fit geometry which repro- 
duces all of the radio peaks, it is clear that a more complex 
radio emission model is needed to explain the observed light 
curve. 

Figure 44 presents the observed and best-fit light curves for 
PSR 12043-^171 1. This is a 2.38 ms pulsar in a binary system 
with a 1.5 day orbital period discovered by Guillemot et al. 
(2012a), who also reported the first detection of gamma-ray 
pulsations from this MSP, in radio observations of unassociated 
LAT sources with pulsar-like characteristics. The OG and TPC 
model both reproduce the observed gamma-ray light curve 
well, though the TPC model matches the peak separation and 
off-peak emission level better. The best-fit geometries using 
both models predict radio profiles that agree qualitatively with 
the observations, but neither can match all of the observed 
components. 

Figure 45 presents the observed and best-fit light curves of 
PSR J2047-F1053. This is a 4.29 ms pulsar in a binary system 
with a 0.12 day orbital period discovered in radio observations 
of unassociated LAT sources with pulsar-like characteristics 
(S. Sanpa-Arsa et al., in preparation). This MSP was first 
announced as a gamma-ray pulsar in 2PC. Both the TPC and OG 
model match the phase of the gamma-ray peak well but neither 
exactly matches the shape. Both models also seem to predict 
the off-peak interval to be from phase 0.8 to 1.1, which does 
not agree well with the data. The best-fit geometries of both fits 
predict two closely spaced radio peaks, though neither exactly 
matches the shape and peak-height ratio. 

Figure 46 presents the observed and best-fit light curves 
of PSR J205 1—0827. This is a 4.51 ms pulsar in a binary 
system with a 0.10 day orbital period discovered by Stappers 
et al. (1996) with eclipses seen for approximately 10% of 
the orbital period at 436 MHz. The observation of eclipses 
suggests that the orbital inclination angle is ^60° . Gamma-ray 
pulsations were first reported from this MSP by Wu et al. (2012) 
and later by Espinoza et al. (2013). Polarimetric observations 
of PSR J205 1—0827 have been reported by several authors 
(Xilouris et al. 1998; Stairs et al. 1999; Ord et al. 2004) who all 
agree that there is some level of linear and circular polarization 
but do not attempt to constrain the viewing geometry. The OG 
and TPC model are both unable to reproduce the gamma-ray 
light curve well, though with current statistics it is unclear 
how many peaks are in the light curve of this MSP. Fits with 
both models do reproduce the observed radio profile well. Only 
the TPC model finds a best-fit ^ > 60°, consistent with the 
observation of radio eclipses. 

Figure 47 presents the observed and best-fit light curves of 
PSR J2124— 3358. This is an isolated 4.93 ms pulsar discovered 
by Bailes et al. (1997). Gamma-ray pulsations from this MSP 
were first reported by Abdo et al. (2009a). Ord et al. (2004), 
Manchester & Han (2004), and Yan et al. (2011) have reported 
polarimetric observations for this pulsar. Manchester & Han 
(2004) attempted RVM fits to the position-angle swing and 
found a = 48° ± 3° and ^ = 67° di 5° but noted a large degree 
of covariance in the parameters, suggesting that a more realistic 
range of values would be ckf g [20°, 60° ] and ^ g [27°, 80°]. 
Yan et al. (2011) argued that if the position angles of linearly 
polarized emission for the peak between —0.3 and —0.2 in 
phase (corresponding to the peak near phase 0.85 in Figure 47) 
were increased by 90°, the position-angle swing would match 
RVM predictions for a ^ 25° and > 0. Bogdanov et al. 
(2008) modeled the thermal X-ray light curve of this MSP 
but the fits were unconstrained in a and The PSPC model 
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is able to reproduce the observed gamma-ray light curve of 
PSR J2124— 3358 well, though the trailing shoulder on the 
main peak is too high in the model. However, our hollow-cone 
beam model is unable to reproduce all the features observed 
in the radio profile, though we do match the peaks at phases 
0.10 and 0.85 well. Our best-fit geometry does predict ^0 > 0 
with a consistent with the prediction of Yan et al. (2011) when 
systematic uncertainties are considered. 

Figure 48 presents the observed and best-fit light curves of 
PSR J22 14-^3000. This is a 3.12 ms pulsar in a binary system 
with a 0.40 day orbital period discovered by Ransom et al. 
(2011), who also reported the first detection of gamma-ray 
pulsations from this MSP, in radio observations of unassociated 
LAT sources with pulsar-like characteristics. PSR J22 14-1-3000 
is a “black-widow” pulsar but no radio eclipses have been 
reported, this suggests that our line of sight is significantly 
above the orbital plane. The PSPC model is able to reproduce 
the two broad peaks observed in the gamma-ray light curve 
well, though the shapes are not exactly correct. The best- 
fit geometry from this fit matches the observed radio profile 
well (both the locations in phase and peak-height ratio). This 
geometry does predict ^ = 86° zb 2°, which is at odds with 
the lack of radio eclipses in a system which is thought to have 
ablated away a significant fraction of the mass of its binary 
companion. 

Figure 49 presents the observed and best-fit light curves of 
PSR 12215-^5135. This is a 2.61 ms pulsar in a binary system 
with a 0.17 day orbital period discovered in radio observations 
of unassociated LAT sources with pulsar-like characteristics 
(J. W. T. Hessels et al., in preparation). This MSP was first 
announced as a gamma-ray pulsar in 2PC. Both the TPC and 
OG models reproduce the observed gamma-ray light curve well, 
and both match the main radio peak but miss the smaller peak 
near phase 0.5. The optical light curve of this MSP has been 
modeled by Breton et al. (2013) and Schroeder & Halpem 
(2014). The first authors found i = 66° zb 16°. The latter 
authors used two different models, finding i — (51.7^.^^ ^)° and 
(51.6t_2^i)°, significantly different from the results of Breton 
et al. (2013). Both the TPC and OG best-fit values of ^ 
agree well with the results of Breton et al. (2013). The OG 
confidence contours can accommodate values of ^ compatible 
with the best-fit i values of Schroeder & Halpem (2014), with 
a near 90°, at the 99% confidence level, but the TPC contours 
cannot. 

Figure 50 presents the observed and best-fit light curves 
of PSR J224 1—5236. This is a 2.19 ms pulsar in a binary 
system with a 0.15 day orbital period discovered by Keith 
et al. (2011), who also reported the first detection of gamma-ray 
pulsations from this MSP, in radio observations of unassociated 
LAT sources with pulsar-like characteristics. These authors also 
reported polarimetric observations for this pulsar but did not 
attempt to constrain the viewing geometry. The OG and TPC 
model both predict two-peaked gamma-ray light curves with 
the correct separation, but neither model properly matches the 
broadness of the peak near phase 0.8. Both models match the 
radio profile well. 

Figure 51 presents the observed and best-fit light curves 
of PSR J2302-I-4442. This is a 5.19 ms pulsar in a binary 
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system with a 51 day orbital period discovered by Cognard 
et al. (201 1), who also reported the first detection of gamma-ray 
pulsations from this MSP, in radio observations of unassociated 
LAT sources with pulsar-like characteristics. Both the TPC and 
OG model match the second gamma-ray peak well but only the 
TPC model predicts the first peak, though not at the correct 
phase. The OG model also has difficulty matching the off-peak 
emission apparent in the top panel of Figure 5 1 . Neither model 
is able to match all of the observed radio components. Johnson 
(2011) also fit the light curves of this MSP with a similar 
likelihood technique and found results compatible with ours. 
They noted that for the geometry corresponding to the TPC 
fit, substantially increasing the emission altitude of the radio 
cone would make the cone broader and more closely match the 
observed profile while leaving the predicted gamma-ray light 
curve unchanged. This would not explain the radio peak near 
phase 0.65. 

APPENDIX B 

SELECT CONEIDENCE CONTOURS 

This appendix presents several example confidence contours 
for each model class. Eor a given MSP, if more than one fit exists, 
then the range of the x axes (^) are matched but the y axes (o') 
are optimized for each fit. The confidence contours show the 
best-fit positions from the tables in Section 7 as stars with 68% 
(light gray, red in the online version), 95% (dark gray, yellow in 
the online version), and 99% (black, green in the online version) 
confidence contours. 

B.l. Contours for Class I MSP s 

Eigures 52 and 55 show the confidence contours for 
PSR J0023-F0923 and 12047-^1053, respectively, using the TPC 
and OG models. These are examples of confidence contours 
that are not simply connected, with two maxima that are al- 
most reflections across the a = line. Eor MSPs with this 
type of confidence contour, changing the radio uncertainty by 
a factor of two can lead to the best fit jumping from one max- 
ima to the other. This leads to systematic uncertainties which 
seem large (changes on the order of 20° to 30°) but that stay 
within the confidence contours, if not the uncertainties quoted in 
Section 7. 

The confidence contours in Eigure 55 are very elliptical and 
tilted with respect to the a and ^ axes. Thus, in addition to not 
being simply connected, the estimated uncertainties in Tables 6 
and 7 overestimate the true size of the contour containing the 
best fit. 

Eigure 53 presents confidence contours for PSR J1024— 07 19. 
These are examples of confidence contours which cover nearly 
the entire allowed parameter space (the comers with large a with 
small ^ and small a with large ^ are excluded because of the 
requirement that the line of sight intersects the radio cone and 
the small a with small ^ region is excluded for the OG model 
due to only having emission above the NCS). The geometry for 
this pulsar is not well constrained due to its faint gamma-ray 
emission. More statistics will improve our constraints. Radio 
polarization measurements and RVM fits would also constrain 
the allowed parameter space. 


45 


The Astrophysical Journal Supplement Series, 213:6 (54pp), 2014 July 


Johnson et al. 


PSR J0023+0923 



Figure 52. Confidence contours in the a-^ plane for PSR J0023+0923 fit with 
the TPC (top) and OG (bottom) models. The star in each panel marks the best- 
fit geometry given in Tables 6 and 7. The light gray (red in the online version) 
contour is 68% confidence, dark gray (yellow in the online version) is 95%, and 
black (green in the online version) is 99%. 

(The complete figure set (28 images) and color version are available in the online 
journal.) 


PSRJ1024-0719 



Figure 53. Confidence contours in the a-^ plane for PSR J1024— 0719 fit with 
the TPC (top) and OG (bottom) models. The star in each panel marks the best- 
fit geometry given in Tables 6 and 7. The light gray (red in the online version) 
contour is 68% confidence, dark gray (yellow in the online version) is 95%, and 
black (green in the online version) is 99%. 

(The complete figure set (28 images) and color version are available in the online 
journal.) 


PSR J 1614-2230 



Figure 54. Confidence contours in the a-^ plane for PSR J1614— 2230 fit with 
the TPC (top) and OG (bottom) models. The star in each panel marks the best- 
ht geometry given in Tables 6 and 7. The light gray (red in the online version) 
contour is 68% conhdence, dark gray (yellow in the online version) is 95%, and 
black (green in the online version) is 99%. 

(The complete hgure set (28 images) and color version are available in the online 
journal.) 
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Figure 55. Confidence contours in the a-^ plane for PSR J2047+1053 ht with 
the TPC (top) and OG (bottom) models. The star in each panel marks the best- 
ht geometry given in Tables 6 and 7. The light gray (red in the online version) 
contour is 68% confidence, dark gray (yellow in the online version) is 95%, and 
black (green in the online version) is 99%. 

(The complete figure set (28 images) and color version are available in the online 
journal.) 
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PSRJ1810+1744 



Figure 56. Confidence contours in the a-^ plane for PSR J1810+1744 fit with 
the alTPC (top), alOG (middle), and laSG (bottom) models. The star in each 
panel marks the best-fit geometry given in Tables 8, 9, and 10. The light gray 
(red in the online version) contour is 68% confidence, dark gray (yellow in the 
online version) is 95%, and black (green in the online version) is 99%. 

(The complete hgure set (6 images) and color version are available in the online 
journal.) 
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Figure 57. Confidence contours in the a-^ plane for PSR J1959+2048 ht with 
the alTPC (top), alOG (middle), and laSG (bottom) models. The star in each 
panel marks the best-ht geometry given in Tables 8, 9, and 10. The light gray 
(red in the online version) contour is 68% confidence, dark gray (yellow in the 
online version) is 95%, and black (green in the online version) is 99%. 

(The complete hgure set (6 images) and color version are available in the online 
journal.) 



Figure 54 presents the confidence contours for PSR J1614— 
2230. These contours, while not simple shapes, are well de- 
scribed by the estimated uncertainties in Tables 6 and 7 with the 
95% confidence level regions mostly concentrated in one area. 

B.2. Contours for Class 11 MSP s 

Figure 56 presents the confidence contours for PSR J1 810-1- 
1744 using the alTPC, alOG, and laSG models. Similar to 
Figures 52 and 55, the alTPC and alOG contours show maxima 
which are not simply connected; there are three regions with 
95% confidence level contours for the alTPC model. The 
laSG contours are connected, but the shape is not simple 
and not easily described by the estimated uncertainties in 
Table 10. 

Figure 57 presents the confidence contours for PSR J1959-f 
2048. These confidence contours are connected and gen- 
erally simple shapes. Both these contours and those for 
PSR J1810-fl744 show a common case where the alTPC and 
alOG best-fit geometries are similar but the laSG geometry is 
quite different, preferring small values of both a and ^ . 

B.3. Contours for Class 111 MSP s 

Figure 58 presents the confidence contours for PSR J1741-f 
1351 fit with the PSPC model. These contours show maxima 
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Figure 58. Confidence contours in the plane for PSR J1741+1351 fit with 
the PSPC model. The star marks the best-ht geometry given in Table 1 1 . The 
light gray (red in the online version) contour is 68% conhdence, dark gray 
(yellow in the online version) is 95%, and black (green in the online version) is 
99%. 

(The complete hgure set (6 images) and color version are available in the online 
journal.) 
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PSR J2214+3000 
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Figure 59. Confidence contours in the a-i; plane for PSR J22 14+3000 fit with 
the PSPC model. The star marks the best-fit geometry given in Table 1 1 . The 
light gray (red in the online version) contour is 68% confidence, dark gray 
(yellow in the online version) is 95%, and black (green in the online version) is 
99%. 

(The complete figure set (6 images) and color version are available in the online 
journal.) 
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which are not simply connected, similar to Figures 52, 55, and 56 
for Class I and II MSPs. These maxima do not show the same 
level of symmetry as in the other cases in which the models have 
emission coming from narrow gaps. 

Figure 59 presents the confidence contours for PSR J2214+ 
3000. This is among the brighter MSPs in gamma rays leading 
to tighter constraints on the emission geometry. This is an 
example where the uncertainties quoted in Table 11 describe 
the confidence region well. 

APPENDIX C 

SIMULATED EMISSION SKYMAPS 

This appendix presents several example simulated emission 
skymaps, corresponding to the best-fit model parameters, for 
each model class. 

C.l. Emission Skymaps for Class I MSPs 

Figure 60 presents the simulated emission skymaps for 
PSR J0030-f0451 corresponding to the best-fit TPC and OG 
models with a hollow-cone radio beam. Both models have high 
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Figure 60. Simulated gamma-ray (top) and radio (bottom) emission skymaps for PSR J0030+0451 fit with the TPC (left, a = 74°) and OG (right, a = 88°) models. 
The white (green in the online version) dashed line in each panel denotes the best-ht ^ . The color scale increases from black to dark gray to light gray to white (black 
to blue to red to yellow to white in the online version). 

(The complete figure set (28 images) and color version are available in the online journal.) 
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Figure 61. Simulated gamma-ray (top) and radio (bottom) emission skymaps for PSR J0613— 0200 fit with the TPC (left, a = 55°) and OG (right, a = 60°) models. 
The white (green in the online version) dashed line in each panel denotes the best- fit ^ . The color scale increases from black to dark gray to light gray to white (black 
to blue to red to yellow to white in the online version). 

(The complete figure set (28 images) and color version are available in the online journal.) 


values of a and ^ . Only in the OG model does the line of sight 
intersect both radio cones leading to an interpulse. Even though 
emission is only collected from below the null charge surface, 
this geometry would allow for off-peak emission for ^ near 0° 
or 1 80° but the predicted radio profile would be much different, 
if the cone was intersected at all. 

Figure 61 presents the simulated emission skymaps for 
PSR J0613— 0200. The radio model for this MSP included both 
a hollow-cone and a core beam, the latter being obvious in 
the lower level emission extending beyond the cone (compare 
to the lower panels of Figure 60, for instance). With more 
moderate values of a, these skymaps demonstrate how no off- 
peak emission is predicted from the OG model. 

Figure 62 presents the simulated emission skymaps for 
PSR J1 024— 0719. This is an example where the best-fit geome- 
tries for the TPC and OG models are the same. From the lower 
panels, this geometry would predict two radio peaks spaced very 
close to each other, but when we reduce the number of bins they 


merge into one peak (see Figure 24). While there does appear 
to be structure in the radio profile shown in 2PC, it is clear that 
our best-fit ^ is not correct; however, there is structure in the 
observed profile which may be matched by a slightly higher ^ 
which can lead to multiple peaks where the emission level does 
not drop all the way to zero inbetween. 

Figure 63 presents the simulated emission skymaps for 
PSR J2017-f0603. This is another case where the TPC and OG 
best-fit geometries are similar, with low a and high ^ . For the 
radio skymap corresponding to the OG fit, lower-right panel, the 
solid angle correction to the bins near o' = 0° leads to very bright 
pixels in this range making the rest of the cone appear artificially 
dimmer compared to the TPC case with slightly larger a. 

C.2. Emission Skymaps for Class II MSP s 

Figure 64 presents the simulated emission skymaps for 
PSR J1823— 3021A corresponding to the best-fit alTPC, alOG, 
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Figure 62. Simulated gamma-ray (top) and radio (bottom) emission skymaps for PSR J1024— 0719 fit with the TPC (left, a = 66°) and OG (right, a = 66°) models. 
The white (green in the online version) dashed line in each panel denotes the best-fit ^ . The color scale increases from black to dark gray to light gray to white (black 
to blue to red to yellow to white in the online version). 

(The complete figure set (28 images) and color version are available in the online journal.) 


and laSG models. The alTPC and alOG skymaps corresponding 
to the radio models, lower panels, are examples of emission 
occurring at high altitudes in the magnetosphere and only over 
a small extent. The laSG skymaps show an example of the 
radio emission region having a smaller emitting gap width 
than the gamma-ray emission. The laSG geometry appears to 
be a large hollow cone, similar to the radio model used for 
Class I and III MSPs but covering a greater portion of the 
sky map. 

Figure 65 presents the simulated emission skymaps for 
PSR J1959-f2048. The alTPC and alOG skymaps show exam- 
ples of the gamma-ray and radio emission regions having the 
same minimum and maximum altitudes. The laSG skymaps for 
this MSP show instances where the emission pattern on the 
sky is more complex than a hollow cone, contrary to the laSG 
skymaps of PSR J1823-3021A. 


C.3. Emission Skymaps for Class III MSPs 

Figures 66 and 67 show the simulated emission skymaps for 
PSR J1744— 1134 and J2214-f3000 corresponding to the best-fit 
PSPC models. Both gamma-ray emission maps, in the top pan- 
els, demonstrate how this emission model can lead to emission 
over most, if not all, of the pulse phase. The radio emission 
skymap for PSR J1744— 1134, the bottom panel of Figure 66, 
shows how just clipping the cone can lead to a sharp radio 
pulse (the dashed line appears to miss the cone but that is just 
an artifact of the color scale, the emission does extend slightly 
beyond the cone visible with this scaling due to the Gaussian 
profile). The radio emission skymap for PSR J2214-f3000, bot- 
tom panel of Figure 67 demonstrates how the line of sight can 
intersect both cones but still give radio peaks with very different 
intensities. 
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Figure 63. Simulated gamma-ray (top) and radio (bottom) emission skymaps for PSR J2017-I-0603 fit with the TPC (left, a = 34°) and OG (right, ol = 23°) models. 
The white (green in the online version) dashed line in each panel denotes the best-ht ^ . The color scale increases from black to dark gray to light gray to white (black 
to blue to red to yellow to white in the online version). 

(The complete figure set (28 images) and color version are available in the online journal.) 



Figure 64. Simulated gamma-ray (top) and radio (bottom) emission skymaps for PSR J1 823— 3021 A ht with the alTPC (left, ot = 46°), alOG (middle, ot = 42°), and 
laSG (right, a = 78°) models. The white (green in the online version) dashed line in each panel denotes the best-ht The color scale increases from black to dark 
gray to light gray to white (black to blue to red to yellow to white in the online version). 

(The complete hgure set (6 images) and color version are available in the online journal.) 
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Figure 65. Simulated gamma-ray (top) and radio (bottom) emission skymaps for PSR J1959+2048 fit with the alTPC (left, a = 56°), alOG (middle, a = 52°), and 
laSG (right, oi = 46°) models. The white (green in the online version) dashed line in each panel denotes the best-ht The color scale increases from black to dark 
gray to light gray to white (black to blue to red to yellow to white in the online version). 

(The complete figure set (6 images) and color version are available in the online journal.) 






Figure 66. Simulated gamma-ray (top) and radio (bottom) emission skymaps 
for PSR J1744— 1134 fit with the PSPC model (o; = 51°). The white (green in 
the online version) dashed line in each panel denotes the best-fit ^ . The color 
scale increases from black to dark gray to light gray to white (black to blue to 
red to yellow to white in the online version). 

(The complete hgure set (6 images) and color version are available in the online 
journal.) 


Figure 67. Simulated gamma-ray (top) and radio (bottom) emission skymaps 
for PSR J2214+3000 fit with the PSPC model (a = 59°). The white (green in 
the online version) dashed line in each panel denotes the best-ht ^ . The color 
scale increases from black to dark gray to light gray to white (black to blue to 
red to yellow to white in the online version). 

(The complete hgure set (6 images) and color version are available in the online 
journal.) 
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